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PREFACE 

Our  long  term  goal  is  to  identify  multiaxial  phenomenological 
constitutive  relations  for  the  finite  strain,  thermomechanical 
behavior  of  unswollen  and  swollen  fiber-reinforced  elastomers. 
Toward  this  end,  we  originally  proposed  to  perform  combined  finite 
extension  and  torsion  tests  on  fairly  large  elastomeric  specimens 
(15  cm  long  x  5  cm  diameter)  using  our  existing  MTS  BIONIX*  load 
frame;  thermal  properties  were  also  to  be  measured  using  standard 
techniques  so  as  to  facilitate  quantification  of  the  thermal 
histories  during  thermoelastic  testing.  Based  on  discussions  with 
ARL  personnel  (soon  after  notification  of  funding) ,  however,  it 
became  apparent  that  it  was  preferred  that  we  develop  a  capability 
to  perform  multiaxial  thermomechanical  tests  on  small  samples  in 
order  to  evaluate  new  candidate  elastomers  as  they  are  developed. 
Hence,  the  scope  of  the  work  was  so  modified.  Because  our  final 
design  and  subsequent  results  will  be  reported  in  refereed  publica¬ 
tions,  the  primary  purpose  of  this  report  is  to  document  our  design 
process,  that  is,  the  theoretical  ideas  and  developments  that 
motivated  the  overall  design,  the  preliminary  tests  that  verified 
some  early  hypotheses,  the  various  design  alternatives  that  were 
considered,  and  details  on  the  final  system  that  was  assembled.  In 
this  way,  this  final  report  will  complement  subsequent  journal 
articles;  together,  they  will  provide  complete  documentation  of 
this  project. 


INTRODUCTION 

It  is  well  known  that  the  properties  of  elastomers  are  highly 
dependent  on  temperature  level,  ranging  from  below  the  glass 
transition  temperature  Tg  to  above  the  melting  temperature  In 
this  research,  we  focus  on  the  behavior  between  Tg  and  Tjj,,  that  is 
in  the  rubbery  range.  A  distinguishing  feature  of  elastomers  within 
T  €  (Tg,Tj„)  is  their  "high  elasticity",  that  is  their  ability  to 
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undergo  large  elastic  deformations.  Not  surprisingly,  there  is  a 
huge  literature  on  the  finite  strain  thermoelastic  behavior  of 
"rubber-like"  materials.  The  earliest  such  paper  is  that  of  Gough 
(1805),  results  of  which  were  latter  confirmed  and  extended  by 
Joule  (1859).  In  particular,  Gough  and  Joule  showed  that  elastomers 
exhibit  "peculiar"  (i.e.,  in  comparison  to  traditional  engineering 
materials  such  as  metals)  behaviors:  for  example,  rubber  warms  when 
stretched  adiabatically,  and  it  shrinks  when  heated.  Although  a 
thermoelastic  theory  was  proposed  by  Kelvin  (1857)  for  elastomers, 
suitable  theoretical  work  had  to  await  a  more  fundamental  under¬ 
standing  of  the  long-chain  polymeric  structure  of  elastomers  (which 
developed  significantly  from  1926  through  1943  due  to  work  by 
Ostwald,  Mark,  Meyer,  Guth,  Flory,  Treloar,  and  others)  and 
advances  in  modern  continuum  mechanics  (which  developed  largely  in 
the  1950s  and  1960s  due  to  work  by  Rivlin,  Noll,  Coleman,  Green, 
Ericksen,  Truesdell,  and  others) .  See  Treloar  (1975)  and  Truesdell 
and  Noll  (1965)  ,  respectively,  for  a  complete  discussion  of  these 
developments.  More  recently,  phenomenological  finite  strain  thermo¬ 
elastic  theories  were  advanced  by  Chadwick  (1974) ,  Chadwick  and 
Creasy  (1984),  Ogden  (1992)  and  others.  These  will  be  discussed 
below.  It  should  be  noted,  however,  that  the  latter  authors  point 
out  that  the  primary  reason  that  well  accepted  descriptors  of  the 
finite  strain  thermoelastic  behavior  of  elastomers  remain  elusive 
is  the  lack  of  complete  multiaxial  experimental  data. 

It  is  similarly  well  known  that  elastomers  are  moderate 
thermal  conductors,  that  is,  reported  values  of  the  thermal 
conductivity  are  three  or  more  orders  of  magnitude  smaller  than 
those  for  metals.  As  with  most  materials,  however,  the  thermal 
conductivity  of  elastomers  varies  with  the  temperature  level  (e.g., 
Dashora,  1994)  .  Schapery  and  Martin  (1974)  and  Burland  et  al. 
(1989)  also  report  some  findings  on  the  thermophysical  properties 
(e.g.,  thermal  diffusivity)  of  particulate  filled  elastomers. 
Nonetheless,  there  is  relatively  little  data  on  the  thermophysical 
properties  of  elastomers  as  a  function  of  multiaxial  finite  strain, 
swelling,  reinforcement  by  polymeric  fibers,  etc. 

The  goal  of  this  work  was  to  develop  a  theoretically-motivated 
test  system  capable  of  providing  much  needed  multiaxial  thermo¬ 
elastic  and  thermophysical  data  on  elastomers  subject  to  a  broad 
range  of  finite  strains,  temperature  levels,  and  environmental 
conditions.  This  technical  report  begins  with  a  review  of  the 
theoretical  motivation,  and  then  docxoments  the  experimental  design. 

THEORETICAL  FRAMEWORK 

General  Concepts 

Coleman  and  Gurtin  (1967)  proposed  a  theoretical  construct  for 
identifying  general  constitutive  functions  that  can  be  used  to 
quantify  the  behavior  of  a  nonlinear ly  elastic,  heat  conducting 
material.  Their  theory  is  briefly  reviewed  here,  albeit  with  some 
modifications  and  notational  changes.  They  suggested  that  a  general 
thermodynamic  process  can  be  described  by  the  following  fields:  the 
current  mass  density  p,  body  force  b,  internal  heat  generation  g, 
Helmholtz  free  energy  i|t,  entropy  r\,  first  Piola-Kirchhof f  stress  P, 
referential  heat  flux  motion  x  of  material  particles  within  the 


body,  and  temperature  T,  all  of  which  can  vary,  in  general,  as  a 
function  of  the  initial  position  X  within  the  body  and  time  t  (note 
that  b,  if,  g  are  defined  per  unit  mass).  Coleman  and  Gurtin 
suggested  that  p,  Jb,  and  g  can  be  found  by  solving  the  balance 
equations  for  mass,  linear  momentum,  and  energy,  which  in 
referential  form  are, 

detF=-^,  f=-|i 

p  aX 

and  J=detF,  F  is  the  deformation  gradient  tensor,  and  p^  the  mass 
density  in  the  reference  configuration, 

'Vg'P+f>Jy=pJt,  where  P=jy^*t  (2) 

and  Vq  =  0(  )fdX,  (  “  )  -  d^{  )/dt^,  (  )"^  denotes  inverse  of  the 
tensor,  t  is  the  Cauchy  (or  true)  stress,  and 

p^6=P’':P-V^-g^+p^gr,  where  (3) 

€  is  the  internal  energy,  q  is  the  spatial  heat  flux  vector,  and 
denotes  transpose  of  the  tensor.  Coleman  and  Gurtin  further 
noted  that  the  free  energy,  entropy,  stress  and  heat  flux  must  be 
determined  from  the  state  variables  (x  and  T,  and  derivatives 
thereof)  via  constitutive  relations.  For  a  local  elastic  theory 
(e.g.,  rate  independent),  appropriate  state  variables  for  a 
homogeneous  material  are  F,  T  and  V^T.  Hence  by  the  principle  of 
equipresence ,  general  (referential)  constitutive  relations  are 

i|r=<IKP,r,v^r) ,  (p,r,v^r) ,  (4)^ 

P=P (P,  T,  V^D  ,  (P,  r,  V^D  (4)2 

where  the  (  "  )  delineates  between  the  function  and  the  value  of 
the  function  (left  hand  sides) .  This  understood,  the  "hat"  notation 
is  dropped  hereafter.  The  balance  of  angular  momentum  and  the 
second  law  of  thermodynamics  (e.g.,  Claus ius-Duhem  inequality) 
provide  restrictions  on  these  constitutive  functions.  In 
referential  form,  these  laws  are 


FP=P^'F^ 

(5) 

-p„(iir+ilf)  +P^:P-Ag^.V^.TiO. 

(6) 

Equation  6  reveals  that  ^  cannot  depend  on  V^T,  and  that  the  four 
constitutive  functions  are  not  all  independent.  Rather, 


dF^ 


^  dT 


(7) 


Hence,  the  stress  and  entropy  must  also  be  independent  of  the 
temperature  gradient,  and  ve  only  need  to  identify  specific  forms 
of  the  free  energy  and  the  heat  flux  in  order  to  have  a  complete 
constitutive  formulation.  Finally,  the  principle  of  material  frame 
indifference  requires  that  i|i  and  depend  on  P  only  through  U,  a 
pure  deformation,  which  renders  the  constitutive  relations 
invariant  to  rigid  body  motion  (recall  the  polar  decomposition 
theorem,  P=R- U,  where  the  orthogonal  tensor  R  represents  rigid  body 
motion)  .  Because  the  right  Cauchy-Green  deformation  tensor  C  equals 


,  however,  i|i  and  can  be  equivalently  written  in  terms  of  C 
(^F^'F).  Hence  we  need  only  identify  specific  functional  forms  of 

i|f=i|r(C,r)  ,  q^=gjC,T,V^T)  ,  (8) 

and  the  associated  values  of  the  material  parameters.  It  is  useful, 
therefore,  to  note  that 

(9) 

is  easier  to  employ  than  equation  7j^.  Moreover,  if  the  material  is 
initially  isotropic  with  respect  to  the  reference  configuration, 
equation  8  reduces  to 

i|r  =i|r  ( D  ,  g„=g„  ( J^,  J J T,  V^r)  ,  ( 10 ) 

where  the  principal  invariants  of  C  are  given  by 

Ic=trC,  2IIc={trO^-tzC^,  JJJj,=detC.  (11) 

Similar  relations  exist  for  transversely-isotropic  and  orthotropic 
symmetries  (Spencer,  1984) ,  which  may  result  from  embedding  fibers 
in  a  homogeneous  elastomeric  matrix.  Nonetheless,  the  goal  of  this 
project  was  to  design  and  construct  an  experimental  system  for 
identifying  forms  of  the  free  energy  and  heat  flux  consistent  with 
equation  8  via  tests  on  small  samples. 

In  experimental  mechanics,  one  always  seeks  a  boundary  value 
problem  that  yields  the  most  easily  and  reliably  interpreted  data. 
For  finite  strains,  experimental  configurations  of  choice  include 
the  in-plane  biaxial  extension  of  thin  sheets  of  material,  combined 
extension  and  torsion  of  solid  cylinders,  inflation  of  membranes, 
and  combined  extension,  inflation  and  torsion  of  thick-walled 
hollow  tubes  (we  have  experience  with  each  of  these  configurations, 
as  applied  to  testing  biological  tissues) .  Of  these  experiments, 
however,  the  simplest  to  interpret  is  the  in-plane  biaxial  test, 
which  is  the  focus  herein. 

Form  of  the  free  energy 

Ogden  (1992)  presented  an  approach  for  inferring  the  form  of 
i|i  from  biaxial  tests,  but  have  we  developed  an  alternate  approach 
which  exploits  the  near  incompressibility  of  elastomers  under 
isothermal  conditions.  Note,  therefore,  that  isothermal  deforma¬ 
tion  induced  volume  changes  in  elastomers  are  typically  0(10“^), 
whereas  temperature- induced  volume  changes  are  of  0(10"'*)  per  °C 
(see  Wood  and  Martin,  1962) .  Hence,  the  latter  will  be  two  orders 
of  magnitude  greater  for  a  100°C  change.  To  exploit  this,  consider 
the  case  of  an  initially  isotropic  thermoelastic  elastomer  that  is 
first  uniformly  heated  from  a  reference  temperature  to  a  T>Tq, 
and  then  stretched  biaxially  (see  Hximphrey,  1996  for  details  of 
this  case  as  well  as  the  case  of  stretch  followed  by  heating) . 
Moreover,  let  the  homogeneous  temperature-induced  and  load-induced 
deformations  (Figure  1)  be  given  by  F-^  -  Si  and  =  diag[A,i, ^2 / ^3]  / 
respectively,  where  A  and  are  stretch  ratios,  detF^’*!,  but 
detF2=l.  Hence,  one  only  needs  to  measure  A,  due  to  changes  in 
temperature,  and  and  A.2»  in-plane  loads  (with  X3  = 

1/ (A.3^A.2)  )  •  Of  course,  the  total  deformation  is  given  by  F=F2'F^  and 
thus  C=F^‘C2‘Fxi  with  C2-F2^'F2-  Now,  since  detC  =  detCidetC2  = 
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detCj^,  equations  22,  9/  and  lOj^  yield  (using  the  chain  rule) 

where  p  is  a  Lagrange  multiplier  that  enforces  the  constraint 
detC2=l/  and  B=F‘ .  Equation  12  is  the  stress-stretch- temperature 
relation.  In-plane,  quasi-static,  biaxial  stretching  in  the  absence 
of  body  forces  results  in  a  homogeneous  plane  stress  field  (tj^3=0, 
i=l,2,3)  in  the  central  region  of  the  specimen,  thus  equilibrium 
(equation  2^)  simply  requires  that  p  is  a  constant  at  each  F2.  The 
non-zero  components  of  the  Cauchy  stress  are,  therefore. 


tij-2pA^  {I1-X3)  (-^^ +A^X2*0^“)  ' 
tj;=2pA“  . 


These  results  represent  2  equations  in  terms  of  2  yet  unknown 
response  functions  (i.e.,  derivatives  of  ii») ,  and  can  be  solved  in 
terms  of  experimentally  "measurable"  quantities  (A,Xj^,tij),  namely 

( J  JJ  III  ,  r)  = - i (14) 

"  2AHXI-XI)  (X^X^)  (X^-X^) 


p-^]lL(r  JJ  jjj  r)= - ± - ( — lii — - — ^ — ).  (14), 

i^TT^  ^  ^  H  .^AA/■^2  ‘i2\  ,a2  i2\  ,a2  i2\ 


2AMx^x^)  ixi-xi)  ai-xi) 


Although  these  response  functions  depend  on  four  state  variables, 
only  two  of  these  (J^  and  JJ^)  change  with  F2,  and  thus  one  can  vary 
these  independently  (at  each  T)  via  traditional  "constant  invariant 
tests"  (Humphrey  et  al.,  1990;  Rivlin  and  Saunders,  1951),  that  is, 
Iq  can  be  maintained  constant  at  multiple  prescribed  values  while 
IIq  is  varied,  and  vice  versa.  When  there  is  no  initial  heating 
from  Tq  (i.e.,  F^^  ^  I) ,  this  protocol  reduces  to  the  now  classical 
case  derived  by  Rivlin  and  Saunders.  Hence,  aside  from  the  same 
caveats  that  exist  in  isothermal  protocols  (e.g.,  sensitivity  of 
results  to  experimental  errors  at  low  strain,  where  and  11^  are 
nearly  equal) ,  equations  14.  ,  allow  one  to  determine  the  form  of 
the  response  ftinctions  at  mul-tiple  temperatures  directly  from  data. 
Note  that  dif/dlllf,  can  be  found  from  trt=0,  for  That  is, 

■=0.  M5^ 


dllc. 


dlllr 


It  appears  we  are  the  first  to  exploit  this  decomposition  of  F,  and 
thus  the  mechanical  incompressibility  constraint. 

One  final  observation  with  regard  to  finding  the  temperature 
dependence  of  f.  Chadwick  (1974),  and  later  Ogden  (1992),  showed 
that  (by  twice  integrating  the  specific  heat  Cp=-r3^i|»/3r^  with 
respect  to  T,  at  a  constant  F)  equation  8^^  can  be  written  as, 

i|r=ilr(C,r„)^+e(C,r^)  (1-^)  *  f\(c,  O  (l-y)dC  (16) 

•‘o  J- O  ‘^'^O  4 

where  i/  -  e-  r]T  and  C  is  a  dummy  variable.  That  is,  a  form  of 
i^(C,r)  can  be  determined  by  knowing  and  e  from  isothermal 


mechanical  tests  at  as  well  as  Cp(C,T).  Equation  16  was  moti¬ 
vated  by  the  desire  to  delineate  between  entropic  and  energetic 
contributions  to  rubber-elasticity,  which  is  essential  in  the 
formulation  of  structurally  based  constitutive  relations  (i.e.,  via 
statistical  mechanics,  see  Treloar,  1975)  .  We  submit,  however,  that 
although  statistical  models  are  useful  in  principle,  the  many 
simplifying  assumptions  currently  required  to  formulate  such 
relations  severely  compromise  their  interpretive  value.  Hence,  our 
focus  Is  on  phenomenological  (ie,  continuum)  relations  that  are 
useful  for  engineering  design  and  analysis  (e.g.,  finite  element 
analyses)  .  For  this  reason,  observe  that  equation  16  can  be  re¬ 
written  as 

i|r=5i  (C,  T^)  +^2  (c,  r^)  r+Sj  (C,  r^,  D  (i7) 

where  are  scalar-valued  functions  that  can  be  determined 

phenomenologically:  and  ^2  can  be  found  from  isothermal 

mechanical  tests  at  T^,  and  ^3  can  be  determined  from  measurements 
of  Cjj  at  multiple  F  and  T,  ^3  will  depend  primarily  on  T,  indeed  it 
provides  the  only  nonlinear  temperature  dependence  to  the  free 
energy.  Equation  17  is  consistent  with  all  current  thermoelastic 
data  (almost  all  of  which  are  uniaxial) ,  which  reveal  that 
elastomers  exhibit  predominantly  linear  stress-temperature 
responses  at  constant  strains  (e.g.,  Anthony  et  al,  1942).  We 
submit  that  equation  17  will  facilitate  the  quantification  and 
evaluation  of  stress-strain-temperature  relations  and  expressions 
for  the  entropy.  In  particular,  we  will  now  be  able  to  compare 
results  from  isothermal  tests  at  multiple  temperatures  (equations 
14,  15)  with  combined  results  from  isothermal  tests  at  and 

measurements  of  Cp(C,T) . 

Form  for  the  heat  flux 

Fourier's  model  for  heat  conduction  is  a  good  descriptor  for 
most  homogeneous  materials  (despite  a  theoretically  unappealing 
infinite  speed  of  heat  diffusion) .  The  common  form  of  this  model  is 
q=-k{T)VT ,  where  k  is  an  isotropic  scalar  thermal  conductivity  and 
this  (spatial)  equation  assumes  that  k  is  independent 
of  strain,  which  is  not  true  in  general.  Hence,  we  consider  a 
generalized  Fourier's  relation, 

go=-jr{c,r)*V„r,  k(c,T)  =6^{c,t)  i+i^{c,T)c+6^{c,T)  (is) 

where  the  form^  for  the  referential  conductivity  tensor  K  follows 
Mueller  (1985) ,  and  the  general  dependence  on  C  can  be  replaced 
with  the  appropriate  invariants  (see  equation  IO2) .  The  importance 
of  this  approach  (which  has  not  been  exploited  experimentally)  is 
that  one  reduces  the  search  for  a  tensor-valued  function  K(C,T)  to 
a  search  for  three  scalar-valued  functions  S^,  which  is  much 
simpler.  Note,  that  K  cannot  be  measured  directly  via  transient 


^  The  Cayley-Hamilton  theorem  for  second  order  tensors  reveals 
that  this  form  is  a  complete  polynomial  expression. 


experiments,  rather  in  such  tests  it  must  be  inferred  from  the 
thermal  diffusivity  aQ(=K/p^Cp)  . 

Measurement  of  thermal  diffusivity 


The  (referential)  thermal  diffusivity  can  be  "measured”  by 
solving  the  energy  equation  (equation  3^)  for  a  tractable  experi¬ 
mental  configuration.  Clearly,  the  energy  equation  is  easier  to 
solve  for  the  case  of  zero  stress-power  (which  can  be  achieved  by 
maintaining  the  in-plane  stretches  constant  in  a  plane  stress 
configuration)  and  zero  internal  heat  generation.  In  this  case, 
Dillon  (1962)  showed  that  equation  3^  becomes, 

where  Cp=-T-^  (19) 


and  thus  equations  18  and  19  yield  a  general  governing  differential 
equation  from  which  the  thermal  diffusivity  can  be  found.  That  is, 
dr/3t=ao:VQ(VpT)  .  Except  for  the  simplest  experiments  (e.g.,  1-D 
heat  conduction) ,  however,  this  equation  is  most  easily  solved 
numerically  (unlike  the  thermoelastic  equilibrium  equation  which 
was  solved  trivially  for  the  biaxial  stretching  configuration) . 

Although  there  are  several  experimental  approaches  for  finding 
Oq,  it  appears  that  the  flash  method  is  preferred  for  materials  of 
moderate  conductivity  that  may  experience  thermal  damage  if 
maintained  at  supra-normal  temperatures  for  long  times.  That  is, 
the  short-lived,  small  temperature  rise  required  for  the  flash 
method  is  less  likely  to  induce  thermal  degradation  of  the 
material;  moreover,  the  flash  method  is  much  better  suited  for 
combination  with  testing  at  multiple  finite  biaxial  strains  because 
it  is  a  non-contact  method.  Briefly,  then,  note  that  in  the  flash 
method  one  subjects  a  thin  specimen  to  a  short  burst  of  radiant 
energy  on  one  face  and  then  measures  the  associated  temperature 
history  on  the  opposite,  non- illuminated  face.  Introduced  by  Parker 
et  al.  (1961)  ,  the  flash  method  was  originally  used  to  measure  the 
out-of-plane  thermal  diffusivity  by  subjecting  the  specimen  to  an 
uniform  illvimination  (Figure  2a)  .  Later,  however,  Maillet  et  al. 
(1990)  and  Lachi  and  Degiovanni  (1991)  showed  that  the  out-of-plane 
and  one  in-plane  component  can  be  measured  simultaneously  if  the 
specimen  is  illuminated  through  a  circular  aperture  (Figure  2b) . 

In  the  case  of  uniform  illumination  (i.e.,  Parker's  original 
method)  ,  one  can  obtain  an  analytic  expression  for  the  out-of-plane 
component  of  the  diffusivity  in  terms  of  experimentally  measurable 
quantities.  Briefly,  the  1-D  energy  equation  for  heat  flow  in  the 
out-of -plane  (i.e.,  3)  direction  is  (from  equation  19) 


dT 

at 


=  (a 


o) 


dXi  ' 


(20) 


where  and  H  is  the  undeformed  thickness.  An  analytic 
expression  for  the  rise  in  the  rear  surface  temperature  T(H,t)  is 
(Carslaw  and  Jaeger,  1959) 
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r(H,  t)  = 
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where  Q  denotes  the  pulse  of  radiant  energy, 
parameters,  <p  and  <■>  of  use  are 


t)  = 


T(H,  t) 


_  ^^(gp)  33^ 


Two  dimensionless 


(22) 


where  is  the  maximum  temperature  seen  on  the  non- illuminated 
face.  Equation  21  can  thus  be  re-written  as, 

e» 

<|>=l+2  (-l)"exp(-n^co) ,  (23) 

n«i 


which  is  a  convergent  series.  Parker  et  al.  suggested  that  this 
equation  be  solved  for  0=0.5,  which  yields  (d=1.38.  Hence,  the  out- 
of-plane  thermal  diffusivity  can  measured  (at  multiple  strains  and 
ambient  temperatures)  via^. 


(®o^  33“ 


1.38^2 


1/2 


(24) 


where  t2^/2  the  time  required  for  the  non-illuminated  surface  to 
reach  half  of  its  maximum  temperature  (Figure  3) .  Although  this  1-D 
analysis  assumes  insulated  edges,  and  that  the  thermal  diffusivity 
is  independent  of  temperature,  these  are  not  severe  restrictions. 
For  example,  the  temperature  rise  is  typically  small  on  the  non- 
illuminated  face,  for  example  from  3  to  5°C.  Although  our  interest 
is  in  the  simultaneous  measurement  of  three  diagonal  components  of 
Oq,  equation  24  provides  a  way  to  (a)  evaluate  our  experimental 
system  and  numerical  methods  for  determining  the  3-D  components  of 
ag(C,T),  and  (b)  determine  a  good  initial  guess  for  the  nonlinear 
regression  methods  that  will  be  needed  to  find  the  3-D  components 
(see  below) . 


Measurement  of  thermal  conductivity  and  specific  heat 

K(C,T)-p^c^{C,T)aQ(C,T)  hence  the  thermal  conductivity  is 
easily  found  from  the  diffusivity,  provided  that  the  original 
density  and  the  specific  heat  are  known:  for  example,  for  the  1-D 
experiment  just  discussed,  the  out-of-plane  thermal  conductivity 
^33  “  (®o)33Po*^F  where  p^Cf^iQ/HT^)  .  Likewise,  recall  (from  equation 
16) ,  that  identification  of  a  form  of  the  free  energy  (and  thus 
stress  and  entropy,  equation  7)  also  requires  that  the  specific 
heat  be  known.  Hence,  quantification  of  c^  as  a  function  of  finite 
strain  and  temperature  is  fundamental  to  a  full  thermomechanical 
constitutive  formulation.  The  specific  heat  c^{I,T)  is  easily 
measured  using  differential  scanning  calorimetry  (DSC) ,  devices  for 
which  are  commercially  available.  DSC  compares  the  energy  required 
to  raise  the  temperature  of  a  small  (unloaded)  specimen  with  the 
energy  required  to  raise  the  temperature  of  a  similarly  sized 
specimen  of  a  specific  heat  standard.  Current  DSC  devices  cannot  be 
used  to  measure  Cp  as  a  function  of  finite  strain,  however,  thus 


^  All  prior  solutions  were  for  the  spatial  diffusivity  a, 
where  a  =  {F-  ol^‘  F^)  fJ ,  and  thus  one  finds  the  same  relation  except 
that  H  is  replaced  with  h  (where  A,3=h/H)  . 


such  data  are  apparently  not  available. 

Consider,  therefore,  that  the  specific  heat,  Cp=d€/dT,  is  a 
scalar  continuum  property  associated  with  either  a  unit  volume  or 
unit  mass.  In  the  absence  of  stress  power,  it  can  be  estimated  to 
be  the  ratio  of  the  change  in  the  total  internal  energy  to  the 
temperature  rise  AT,  that  is  c^=A€^/jnAr  where  m  is  the  mass  of  the 
sample.  Yet,  this  equation,  which  forms  the  basis  of  DSC  analysis, 
assumes  a  uniform  temperature  change  throughout  the  specimen.  This 
will  not  be  true  in  general.  Indeed,  we  wish  to  explore  the  possi¬ 
bility  of  determining  the  specific  heat  of  a  specimen,  at  various 
states  of  finite  strain  and  temperature  level,  from  a  thin  biaxial 
sample  that  is  subjected  to  a  nearly  uniform  burst  of  radiant 
energy.  In  such  a  case  the  temperature  in  the  specimen  will  clearly 
be  non-uniform  in  the  in-plane  dimensions;  the  out  of  plane 
dimension,  being  at  most  1/30^*^  of  those  in-plane,  will  reach 
thermal  equilibrium  more  quickly.  Hence,  consider  a  volumetrically 
averaged  specific  heat, 

ATjA  V,)  -1  (25) 

where  "pixels”  suggests  that  the  temperature  field  on  the  entire 
non- illuminated  face  will  be  known  (e.g.,  digitized  from  the  output 
of  an  infrared  camera)  .  At  present,  however,  this  is  only  a 
preliminary  idea;  much  work  remains  to  determine  the  best  way  to 
determine  c^{C,T)  . 

A  FEASIBILITY  EXPERIMENT 

In  conjunction  with  our  theoretical  developments,  we  assembled 
a  simple  biaxial  loading  frame  to  test  our  hypothesis  that  the 
dependence  of  on  C  cannot  be  neglected  in  the  elastomers  of 
interest  here.  Constructed  of  1.9  cm  thick  303  stainless  steel,  the 
17.8  X  17.8  cm  loading  frame  (Figure  4)  could  be  placed  within  a 
30.5  cm  diameter  bell  jar  so  that  experiments  could  be  performed 
under  vacuxim  to  evaluate  possible  effects  of  convection.  Oil- 
impregnated  bronze  bushings  guided  orthogonally  mounted  pairs  of 
stainless  steel  loading  carriages,  which  in  turn  subjected  square 
specimens  to  finite  biaxial  extensions  via  four  arrays  of  loading 
strings.  Specimens  were  loaded  by  hand  by  turning  four  bolts,  and 
allowed  to  stress  relax  for  10  minutes.  Associated  extensions  were 
measured  using  a  hand  caliper,  but  the  loads  were  not  measured  — 
we  were  only  interested  in  determining  the  effects  of  finite  strain 
on  Oq,  not  in  constitutive  formulations. 

A  flash  system  was  developed  (described  below)  and  used  to 
uniformly  illuminate  (via  typically  50  J  in  135  /is)  multiple 
silicone  samples.  Initial  measurements  of  T  on  the  non-illuminated 
face  were  accomplished  using  a  "thermocouple  probe".  This  probe 
consisted  of  two  0.254  mm  diameter  T-type  (copper-constantan) 
thermocouples,  each  of  which  were  inserted  into  a  0.794  mm  diameter 
ceramic  insulator.  The  probe  was  then  mounted  on  a  1/4"  x  20 
stainless  steel  machine  screw  which  allowed  the  thermocouples  to  be 
brought  into  contact  with  the  specimen  (with  an  unknown  contact 
load)  by  advancing  a  screw.  These  small  diameter  thermocouples  had 
a  rapid  response  time  (0(0.05)  s)  in  still  air,  but  we  switched  to 
0.51  mm  diameter  E-type  (chromel-constantan)  thermocouples  which 


have  the  highest  Seebeck  coefficient  in  the  range  of  temperatures 
considered.  This  E-type  probe  was  used  to  gather  most  of  the 
initial  data,  and  is  still  being  used  (albeit  modified  to  hold 
three  thermocouples) .  Temperatures  were  initially  monitored  using 
a  Fluke  Hydra  model  2620-A;  this  unit  is  only  capable  of  monitoring 
three  thermocouples  every  1.3  seconds,  with  an  approximate  uncer¬ 
tainty  of  0.5°C.  For  1.5  mm  thick  silicone  specimens,  the  t2/2 
(equation  24)  was  about  4  seconds  and  was  about  3°C.  Hence,  this 
system  did  not  have  adequate  temperature  or  time  resolution.  Next, 
an  available  Data  Translation  A/D  board  (DT  2805)  was  used  (in  a 
386SX  IBM  compatible  computer)  to  try  to  get  better  resolution. 
This  system  was  capable  of  measuring  three  channels  in  0.143 
seconds,  but  only  to  within  approximately  0.5°C.  This  system  was, 
therefore,  also  inadequate.  Nonetheless,  associated  pilot  data 
collected  in  March  1995  revealed,  as  expected,  that  the  out-of¬ 
plane  diffusivity  (calculated  using  equation  24)  varied  signifi¬ 
cantly  with  finite  strain.  Figure  5  shows  a  representative 
temperature  rise  on  the  non- illuminated  face,  and  Figure  6  shows 
results  for  the  out-of-plane  diffusivity  as  a  function  of  finite 
uniaxial  strain.  These  pilot  data,  combined  with  aforementioned 
theoretical  considerations,  justified  the  design  of  a  new  system 
capable  of  determining  multiaxial  thermoelastic  and  thermophysical 
properties  over  a  broad  range  of  finite  strains,  temperatures,  and 
environmental  conditions.  This  system  is  described  in  detail  below. 

NUMERICAL  SIMULATIONS 

Finite  element  stress  analysis 

Equations  13-15  and  20  assume  that  the  stress  and  strain 
fields  are  homogeneous,  that  is  independent  of  position  within  the 
sample.  In  an  actual  experiment,  however,  there  will  be  significant 
variations  in  both  fields  because  of  the  necessity  to  load  the 
samples  with  multiple  point  loads  as  seen  in  Figure  4  (this  allows 
the  specimen  to  deform  freely  in  the  direction  orthogonal  to  the 
applied  load) .  Though  it  has  been  shown  empirically  that  strains 
are  homogeneous  within  a  central  region  of  biaxial ly  loaded  samples 
(Treloar,  1948;  Humphrey  et  al.,  1990) ,  and  via  finite  elements  for 
one  class  of  elastomers  (Nielson  et  al.,  1991),  we  performed  a 
finite  element  analysis  to  determine  the  extent  of  the  central 
homogeneous  field.  Assuming  a  Mooney-Rivlin  type  hyperelastic 
behavior  (i.e.,  W  =  Ci{Iq-3)  +  C2(I’Jc“3),  where  are  material 
constants  and  W  is  the  strain  energy) ,  the  ABAQUS®  finite  element 
code  was  used  to  study  possible  experimental  configurations,  that 
is  the  effect  of  the  position  and  nxamber  of  point  loads.  Because  of 
double-symmetry,  one-quarter  of  the  specimen  was  studied;  1200 
(i.e.,  20x20x3)  hybrid  brick  elements  were  used  to  discretize  the 
domain,  which  was  taken  to  be  15  mm  per  side.  Cases  were  run  with 
even  or  odd  numbers  (4, 5, 6, 7  and  8)  of  equally  spaced  point  loads 
per  side.  Figures  7a  and  7b  show  results  for  an  equibiaxial  stretch 
of  10%.  As  expected,  we  found  that  at  least  5  point  loads  were 
needed  per  side  to  ensure  that  the  stress  and  strain  fields  within 
the  central  25%  of  the  overall  area  was  homogeneous  to  within  3%. 

Thermal  diffusivity  calculations 


As  mentioned  above,  determination  of  the  three  (diagonal) 
components  of  requires  that  the  3-D  heat  equation  be  solved  for 
the  temperature  field  T(X,t)  throughout  a  biaxially  stretched 
specimen  that  is  exposed  to  a  burst  of  radiant  energy  over  part  of 
the  "illxuninated  surface”  (recall  Figure  2b) ;  indeed,  multiple 
solutions  are  needed  at  various  fixed  C  and  various,  constant 
ambient  temperatures.  Clearly,  this  solution  must  be  achieved 
numerically,  the  accuracy  of  which  depends  primarily  upon  how  well 
boundary  conditions  are  approximated.  It  is  usually  assumed  in  the 
flash  method  that  a  uniform  and  instantaneous  heat  pulse  is 
absorbed  in  a  thin  layer  of  a  homogeneous,  opaque,  and  insulated 
sample.  In  reality  these  conditions  are  seldom  met,  especially  for 
samples  of  low  to  moderate  thermal  conductivity.  In  order  to 
understand  these  problems,  2-D  models  have  been  used  in  which 
thermal  diffusivity  has  been  determined  by  temporal  moments 
(Digiovanni  and  Laurent,  1986)  or  by  use  of  Laplace  transforms 
(Gembarovic  and  Taylor,  1993).  Nonlinear  regression  schemes  have 
also  been  used  extensively  (Gembarovic  et  al.,  1990;  Vozar  et  al. , 
1991;  Sawaf  and  Ozisik,  1995) .  For  example,  Sramkova  and  Log  (1995) 
used  the  Marquardt-Levenberg  estimation  scheme  in  combination  with 
a  model  that  accounted  for  heat  loss.  They  provided  a  comparison  of 
results  from  other  authors  techniques  using  NIST  standard 
materials,  and  claimed  excellent  results. 

Our  approach  has  been  to  combine  the  Marquardt  regression 
method  with  a  finite  difference  solution  of  the  heat  equation.  In 
particular,  the  heat  equation  is  solved  using  central-differencing 
in  space  and  explicit  first  order  differencing  in  time.  An  explicit 
method  offers  several  advantages  over  an  implicit  method.  For 
example,  the  implicit  model  requires  simultaneous  solution  of  all 
nodal  temperatures  at  each  time  step  —  this  requires  inverting  a 
large  matrix  at  each  time,  thus  this  method  is  time  consuming  when 
the  number  of  discrete  time  points  is  large.  In  contrast,  an 
explicit  method  calculates  present  time  directly  from  past  time; 
because  the  solution  marches  forward  in  time,  it  does  not  require 
as  much  computational  time  to  calculate  temperatures  at  each 
discrete  time.  One-,  two-,  and  three-dimensional  codes  have  been 
written  in  FORTRAN  to  analyze  modified  flash  experiments.  Before 
discussing  results  from  our  preliminary  simulations,  however,  first 
consider  the  parameter  estimation  method. 

The  Marquardt-Levenberg  algorithm  has  been  shown  to  be  robust 
for  determining  elastic  parameters  from  biaxial  stretching  tests  on 
rubber  (Twizell  and  Ogden,  1983)  and  tissues  (Humphrey  et  al., 
1990) ,  as  well  as  thermal  conductivities  in  solids  (Sawaf  and 
Ozisik,  1995) .  Briefly,  this  algorithm  is  a  nonlinear  least-squares 
method  wherein  one  seeks  to  minimize  the  objective  function  e, 

®=E?=i  (yt(P)-ye)i  (26) 

where  Jij^  is  the  total  number  of  measurements  (i.e.,  number  of 
measurement  sites  times  the  number  of  data  points  at  each  site) , 
the  vector  of  theoretically  determined  values  (e.g.,  stresses  or 
temperatures),  p  the  vector  of  material  parameters  (e.g.,  material 
moduli  or  dif fusivities) ,  and  the  vector  of  experimentally 
determined  values.  Values  of  p  are  found  iteratively,  via 


(27) 


where  i  is  an  iteration  counter  and  v  the  Marquardt  parameter;  v  is 
typically  initialized  to  “0.01,  and  multiplied  or  divided  by  10  if 
e  increased  or  decreased  during  the  prior  iteration..  Iterations 
continue  until  a  norm  of  the  difference  between  the  and  p^^^ 
vectors  becomes  less  than  a  prescribed  value.  We  have  extensive 
experience  with  multivariate,  constrained  regressions  using  the 
Marquardt  method  (e.g.,  see  Humphrey  et  al.,  1990). 

To  determine  the  three  diagonal  components  of  the  thermal 
diffusivity,  we  must  minimize  the  difference  between  the  tempera¬ 
tures  measured  and  calculated  at  three  different  locations  on  the 
non-illuminated  face  of  a  specimen  (see  Figures  8a  and  8b)  at 
multiple  times  (i.e.,  not  just  at  in  the  1-D  problem, 
equation  24)  .  Furthermore,  the  Marquardt  regression  must  be 
initialized  by  providing  initial  guesses  for  the  values  of  each  of 
the  parameters,  that  is  p^^^;  it  is  well  known  that  the  Marquardt 
method  is  relatively  insensitive  to  the  initial  guess,  but 
certainly  the  better  the  guess  the  more  reliable  and  faster  the 
regression.  In  practice,  good  initial  guesses  can  be  obtained  by 
using  equation  24  and  a  uniform  flash  to  determine  the  out-of -plane 
component,  and  then  using  this  value  as  a  guess  for  all  three 
components  in  a  full  regression,  given  the  same  biaxial  strain 
state  and  ambient  temperature.  Best-fit  values  for  each  state  can 
then  be  used  as  initial  guesses  for  subsequent  states  (i.e.,  strain 
and  temperature  level) . 

Before  evaluating  the  combined  Marquardt-Finite  Difference 
code,  however,  we  first  verified  that  the  finite  difference 
solution  was  reasonable.  For  example,  we  used  a  2-D  model  to 
evaluate  the  uniform  flash  experiment  for  which  a  1-D  solution 
(equation  24)  is  available.  The  "specimen"  was  taken  to  be  10  x  10 
X  0.1  cm,  and  assxamed  to  be  silicone  rubber  (HSII)  having  an 
isotropic  thermal  conductivity  a  of  0.22  W/m^K.  The  specimen  was 
exposed  to  50  J  of  energy  with  a  flash  time  of  100  ns.  Figure  9 
shows  the  non-illuminated  face  temperature  as  a  function  of  the 
position  from  the  centerline  of  the  sample;  temperature  was  non- 
dimensionalized  as 


AT 


(28) 


where  is  the  absorbed  energy,  V  the  specimen  volume,  and  the 
overbar  denotes  a  volumetric  average.  As  can  be  seen  from  Figure  9, 
>95%  of  the  non-illuminated  face  of  the  sample  achieved  a  uniform 
temperature  (after  4t3^y2)  that  was  >90%  of  T^.  That  is,  losses  at 
the  edge  of  the  sample  were  not  significant  within  the  observation 
period  of  interest,  and  therefore  the  1-D  solution  is  appropriate. 
A  1-D  explicit  finite-difference  model  for  the  uniform  flash 
method,  including  convection  effects,  is  currently  being  compared 
to  the  analytic  1-D  solution  for  further  validation  of  the 
numerical  solution.  Nonetheless,  the  importance  of  numerical 
simulations  is  evident;  this  study  revealed  that  uniform  flash 
experiments  are  tractable  for  homogeneous  rectangular  elastomeric 


specimens,  and  convection  losses  are  minimal  even  when  the 
experiments  are  performed  in  air  rather  than  in  a  vacuum. 

A  3-D  explicit  Cartesian  finite-difference  model  was  then 
developed  to  study  the  non-uniform  illumination  experiment:  that 
is,  illumination  of  a  square  specimen  through  a  smaller  square 
aperture  as  seen  in  Figure  8a  (recall  from  Figure  2b  that  Malliet 
et  al.  (1990)  used  a  circular  aperture  and  a  square  specimen) . 
Figure  10  reveals  the  model  and  boundary  conditions:  the  aperture 
size  was  6.5  x  6.5  mm  and  the  specimen  size  was  30  x  30  mm,  but 
symmetry  allowed  examination  of  only  one  quadrant.  The  boundary 
conditions  were  (a)  insulated  centerline  edges,  (b)  natural 
convection  on  the  outer  edges  and  non-illuminated  face,  and  (c)  a 
step  from  a  constant  heat  flux  to  a  convective  boundary  on  the 
illuminated  face.  Figure  11a  shows  representative  temperature 
histories  at  different  positions  on  the  X2  axis  of  the  specimen. 
This  model  revealed  locations  for  thermocouples  on  actual  specimens 
that  will  yield  a  measurable  temperature  rise  in  the  in-plane 
directions.  Figure  lib  reveals  further,  however,  that  the  tempera¬ 
tures  calculated  in  the  in-plane  directions  exhibited  a  strong 
nodal  dependence,  likely  due  to  steep  temperature  gradients  at  the 
boundary  of  the  illiaminated  region.  For  example,  there  was  a  3.7°C 
change  in  the  in-plane  directions  in  comparison  to  a  0 . 1°C  change 
in  the  out-of-plane  direction  over  the  same  difference  in  nodes 
(i.e.,  10  and  100  nodes).  Hence,  the  mesh  must  be  chosen  carefully 
to  ensure  reliable  parameter  estimation. 

Given  that  the  finite  difference  model  gave  reasonable 
results,  we  then  performed  a  sensitivity  study  to  evaluate  the 
robustness  of  the  combined  Marquardt-Finite  Difference  code. 
Briefly,  we  prescribed  values  of  the  diffusivity  and  specimen 
dimensions,  and  used  the  finite  difference  algorithm  to  solve  the 
temperature  field  associated  with  a  prescribed  "flash”.  Values  at 
three  locations  on  the  non-illuminated  face,  that  is  T{X^,X^,H ,t) , 
were  then  perturbed  by  adding  random  Gaussian  noise  (to  simulate 
experimental  errors)  and  used  as  input  to  our  Marquardt  regression 
as  experimental  data.  In  this  way,  we  could  evaluate  the  robustness 
of  our  algorithm  as  a  function  of  the  location  of  the  temperature 
measurements,  number  of  data  points,  amount  of  experimental  noise, 
goodness  of  our  initial  guess,  etc.  As  expected,  the  combined 
Marquardt-Finite  Difference  parameter  estimation  worked  well  for 
the  1-D  model  (for  which  needs  to  determine  only  one  parameter)  . 
Nonetheless,  we  are  now  comparing  the  finite  difference  and 
analytical  based  (equation  24)  regressions  for  speed  and  accuracy. 
The  results  from  these  tests  will  further  support  the  use  of  the 
niimerical  model  in  our  experimental  data  reduction.  Likewise,  the 
3-D  explicit  finite  difference  model  was  incorporated  into  a 
multivariate  Marquardt  parameter  estimation  scheme  to  deteirmine  the 
three  diagonal  components  of  a^.  Initial  evaluation  of  this  code, 
using  numerically  generated  "data"  without  noise,  revealed  that  we 
could  recover  known  parameters  to  within  0.1%  in  less  than  10 
iterations.  As  expected,  convergence  was  quicker  in  the  out-of- 
plane  direction  than  in  the  in-plane  directions.  This  is  due  to 
less  overall  sensitivity  of  the  thermal  diffusivity  in  the  out-of- 
plane  direction,  thus  thermocouple  location  and  aperture  alignment 
have  a  substantial  effect  on  the  determination  of  in-plane  thermal 
diffusivities.  Additional,  more  complete  estimations  are  currently 


in  progress. 

EXPERIMENTAL  SYSTEM 

Figure  12  is  a  schema  of  the  overall  optical-thermomechanical 
system.  Briefly,  the  system  consists  of  six  sub-systems,  that  is 
one  each  for  (a)  biaxial  loading,  (b)  in-plane  biaxial  strain 
measurement,  (c)  thickness  measurement,  (d)  ambient  temperature 
control,  (e)  flash  illumination,  and  (f)  point-wise  temperature 
measurement.  Each  sub-system  is  discussed  in  order  below. 

Biaxial  loading 

Oblique  and  top  views  of  the  device  are  shown  in  Figures  13a 
and  13b.  The  loading  frame  (LF)  is  constructed  from  one  piece  of 
mild  steel  having  overall  dimensions  of  45.72  x  45.72  x  2.54  cm, 
and  a  2.54  x  2.54  cm  cross  section.  The  frame  is  mounted,  at  each 
corner,  on  an  optical  table  using  four  standard  2.54  cm  diameter, 
15.24  cm  long  optical  rods  to  ensure  stability  (Figure  13a).  All 
surfaces  in  the  plane  of  the  frame  are  ground  so  as  to  ensure  flat 
and  parallel  faces  for  mounting  the  various  components  that  make  up 
the  loading  frame.  Pairs  of  through  holes  (1.27  cm  diameter  and 
4.445  cm  apart)  house  recirculating  bearings  (1.905  cm  long)  on 
each  side  of  the  frame.  These  bearings  allow  linear  travel,  with 
minimal  play,  of  four  pairs  of  0.635  cm  diameter  solid  steel 
loading  shafts,  which  form  two  opposing  sets  of  load  carriages 
(CAR) .  Each  shaft  is  also  supported  by  a  rulon  bearing  that  forms 
a  watertight  seal  within  the  wall  of  the  environmental  chamber 
(EC) .  Proper  alignment  of  the  environmental  chamber  and  loading 
frame  ensure  minimal  sticking  and  thus  smooth  translation  of  the 
load  carriages. 

The  inside  half  of  the  load  carriage  consists  of  four  separate 
parts  (Figure  14)  ,  all  made  of  stainless  steel.  The  first  part  is 
a  cross  bar  that  serves  to  hold  the  two  shafts  in  place  as  well  as 
holding  a  load  cell  (LC)  or  a  spacer  rod.  The  load  cells  are  easily 
interchangeable,  but  are  currently  9.8  N  Sensotec  transducers  that 
are  temperature  compensated  (over  the  range  -2 OF  to  200  F)  and 
shielded  against  wet  environments.  These  cells  will  be  used  during 
tests  on  natural  rubber,  whereas  445  N  load  cells  will  be  used  for 
tests  on  other  polymers,  such  as  polyurethane.  The  load  cells  are 
excited  by  a  signal  conditioner /amplifier  (SA-4)  ,  which  in  turn  is 
connected  to  a  12  bit  A/D  board  for  load  registration.  Mechanically 
attached  to  the  load  cells  and  spacer  rods  is  a  coupling  bar  which 
holds  an  interchangeable  T-section.  This  T  has  a  series  of  equally 
spaced  holes  (e.g.,  7.5  mm  apart)  which  will  serve  as  rigging 
points  for  the  arrays  of  Kevlar*  (MIL-T-87128)  threads  (K)  that 
attach  to  the  specimen  (S) .  On  the  other  end  of  the  shafts  is  a 
similar  cross  bar  that  holds  the  shafts  in  a  similar  fashion,  and 
facilitates  load  application. 

Currently,  the  load  is  applied  by  hand  using  a  threaded  rod. 
That  is,  a  15.24  cm  long  1/2 ”-13  threaded  rod  is  attached  to  a  face 
plate  that  is  mounted  directly  on  the  outer  surface  of  the  load 
frame.  The  rod  passes  through  a  1.27  diameter  hole  in  the  outside 
crossbar.  Linear  translation  of  the  load  carriage  is  achieved  by 
manually  rotating  a  nut  on  the  inside  of  the  crossbar  in  either 


direction.  Although  this  system  is  sufficient  for  many  of  our 
envisioned  tests  (which  require  both  thermal  and  mechanical  equili¬ 
brium)  ,  we  are  in  the  process  of  converting  to  micro-stepper  motor 
control.  Four  independently  controlled  motors  (M)  will  drive  the 
loading  rods  via  a  ball  screw  and  ball  nut  (attached  to  the  cross 
bar)  .  The  motors  will  be  attached  to  the  ball  screw  using  flexible 
couplings,  and  supported  from  the  tabletop  via  a  bracket  (SUP) .  By 
using  two  pairs  of  oppositely  moving  loading  frames,  the  center  of 
the  specimen  will  remain  fixed,  which  in  turn  facilitates  strain 
measurement.  The  four  motors  will  be  controlled  individually  via  a 
four  axis  indexer  and  controller  residing  in  a  PC;  we  have 
considerable  experience  with  similar  devices  (e.g.,  see  Humphrey  et 
al.,  1993). 

Biaxial  strain  measurement 

Because  of  the  homogeneity  in  the  central  region  of  the 
specimen,  in-plane  strains  (extension  and  shear)  are  easily 
calculated  from  the  motions  of  3  to  4  surface  markers.  We  have  used 
this  technique  extensively  before,  which  is  described  elsewhere 
(Humphrey  et  al.,  1990;  Downs  et  al.,  1990).  Briefly,  however,  note 
that  four  small  white  markers  will  be  affixed  to  the  central  region 
of  the  specimen  on  the  face  that  will  be  illuminated  during  a  flash 
(recall  Figure  8b)  —  the  specimen  will  be  darkened  in  order  to 
improve  the  absorption  of  radiant  energy,  thus  white  markers  will 
contrast  well  against  the  dark  background  and  thereby  be  easily 
tracked.  Note,  too,  that  strains  will  be  measured  continuously 
during  thermoelastic  tests  (in  the  absence  of  a  flash)  and  at 
mechanical  and  thermal  equilibrium  prior  to  a  flash.  Hence  there 
will  be  no  interference  between  the  flash  and  video-based  strain 
measurement.  Having  the  strain  and  flash  systems  focus  on  the  same 
face  also  leaves  the  non- illuminated  face  (actually  the  space  above 
it)  free  for  the  thermocouple  probe.  Nonetheless,  the  centroids  of 
the  markers  will  be  tracked  on-line  at  the  30  Hz  frame  rate  using 
custom  correlation-based  software  (Downs  et  al.,  1990).  Given  the 
marker  positions,  one  can  use  interpolation  functions  (e.g., 
bilinear)  to  construct  continuous  displacement  fields  within  the 
central  region,  which  in  turn  can  be  used  to  calculate  the 
displacement  gradient  tensor  H,  Since  F=J+H,  the  deformation 
gradient  can  thereby  be  determined,  and  hence  all  necessary 
measures  of  strain. 

The  marker  motions  will  be  monitored  with  a  CCD-Iris  camera 
(Sony  SSC-M370)  via  an  optical  mirror  (mounted  directly  beneath  the 
specimen,  and  at  45°  to  the  optical  table) ;  the  CCD  and  reflecting 
mirror  are  denoted  by  C  and  R,  respectively,  in  Figure  13b.  The  CCD 
has  a  1.27  cm^  chip  area,  and  its  output  is  fed  directly  into  a  DT 
2853  video  frame  grabber  board  as  well  as  a  B&W  (Panasonic  TR-990C) 
monitor  (Figure  12) .  The  frame  grabber  digitizes  the  image  in  real 
time,  thereby  allowing  the  correlation  algorithm  to  determine  the 
marker  locations.  Currently,  we  are  using  a  Canon  100  mm  f4  macro 
lens,  which  has  a  magnification  of  0.5  at  its  focal  length.  A  C- 
mount  adaptor  couples  the  lens  to  the  CCD. 


Thickness  measurement 


Although  the  mechanically-induced  deformations  are  nearly 
isochoric  (i.e.,  A,3=l/>.j^A.2/  see  equation  12  and  the  related 
discussion) ,  the  thermally  induced  deformations  are  not  (recall 
Figure  1  and  that  F’3=AJ)  .  Hence,  the  current  thickness  h  of  the 
specimen  must  be  measured  as  temperature  changes.  We  have 
investigated  several  options  for  on-line  thickness  measurement,  but 
we  have  yet  to  find  a  satisfactory  method.  One  option,  however, 
which  is  the  simplest  in  principle,  is  to  recall  that  detf^A^^p./p, 
where  the  current  mass  density  p  will  change  with  temperature  — 
volume  changes  are  on  the  order  of  per  degree  C.  Hence,  by 
simply  measuring  p(r),  the  temperature-induced  thickness  changes 
can  be  calculated.  This  is  being  explored. 

Environmental  chamber  and  ambient  temperature  control 

As  mentioned  above,  we  seek  to  measure  thermoelastic  and 
thermophysical  properties  at  (ambient)  temperatures  levels  over  the 
range  T  e  (20,100°C],  and  under  a  variety  of  environmental 
conditions  (e.g.,  humidity,  solvents,  etc.).  This  requires  a  closed 
environmental  chamber  and  active  temperature  control.  The  chamber 
(EC  in  Figure  13)  is  a  layered  sandwich  structure  constructed  of  a 
thin  inner  layer  of  silicone  rubber  (to  prevent  leakage)  ,  a  central 
layer  of  polycarbonate  (for  stiffness) ,  and  an  outer  layer  of 
polystyrene  for  insulation  against  parasitic  heat  losses  due  to 
conduction  and  radiation.  Specifically,  the  1.27  cm  thick  central 
polycarbonate  shell  is  40  x  40  x  7.62  cm  in  outer  dimension,  and 
thus  fits  snugly  on  the  inside  of  the  load  frame.  A  pair  of  holes 
on  each  side  of  the  chamber  allow  the  loading  rods  to  pass  through 
to  the  inside.  These  holes  are  fitted  with  neoprene  rubber  seals, 
consisting  of  a  rulon  bearing  and  an  o-ring  at  each  end  of  the 
bearing;  this  minimizes  any  leakage,  and  reduces  heat  loss  (except 
due  to  conduction  along  the  rods) . 

Two  face  plates,  one  on  the  top  and  one  on  the  bottom,  are 
attached  to  the  shell  using  8  nylon  bolts  each;  0.16  cm  thick 
silicone  gaskets  prevent  heat  loss  and  fluid  leakage  from  the 
chamber.  Note,  too,  that  each  face  plate  is  designed  differently. 
The  bottom  plate  (Figure  13a)  contains  a  centrally  mounted,  thin 
quartz  window  through  which  the  specimen  can  be  heated  using  the 
flash  system  (see  below)  or  monitored  using  the  video-based  strain 
measurement  sub-system.  The  quartz  window  is  raised  relative  to 
majority  of  the  surface  of  this  face  plate,  and  thus  is  1.27  to 
1.905  cm  from  the  bottom  surface  of  the  specimen.  This  allows  the 
flash  lamp  to  be  placed  closer  to  the  specimen,  which  in  turn 
improves  heat  transfer.  The  bottom  face  plate  also  has  two  drainage 
ports  to  remove  fluid  from  the  chamber  when  needed  -  recall  that 
some  thermoelastic  tests  will  be  performed  with  the  specimen 
immersed.  Adequate  space  is  provided  under  this  face  plate  for  both 
an  optical  mirror  (for  the  strain  measurement  system)  that  is 
mounted  on  the  optical  table,  and  the  flash  system  which  is  mounted 
to  the  plate  (F  in  Figure  14) .  The  top  face  plate  is  much  simpler. 
Besides  the  necessary  access  ports  for  wiring  etc.,  the  only  other 
feature  is  an  2.54  cm  diameter  central  hole  which  allows  the 
thermocouple  probe  to  contact  the  central  region  of  the  non- 
illuminated  face  at  select  locations  on  the  top  surface  of  the 
specimen  (recall  Figure  8) . 
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Prior  to  the  final  design,  the  heat  loss  from  the  chamber  was 
estimated  for  the  maximum  expected  ambient  T  =  100°C.  This  was  done 
by  taking  into  account  the  quartz  windows  and  access  port  for  the 
thermocouple  probe  as  well  as  conduction  and  radiation  losses 
through  the  walls  and  the  support  rods  of  the  carriages.  The  total 
loss  was  calculated  (using  1-D  resistance  models  for  the  walls)  to 
be  approximately  700  W.  As  expected,  the  estimated  heat  loss  from 
the  walls  of  the  chamber  and  the  support  rods  of  the  load  carriages 
was  1/10  that  through  the  windows  which  are  not  covered  with  the 
2.54  cm  thick  polystyrene.  Hence,  to  maintain  ambient  equilibrium 
temperatures  up  to  100°C,  and  some  symmetry  in  the  heating,  power 
is  to  be  supplied  by  a  pair  of  immersible  rod  heaters  (H)  that  fit 
within  the  confines  of  the  chamber  (Figure  13) .  These  rods  (Omega) 
are  1.27  cm  in  diameter  and  30.5  cm  long,  and  require  a  maximum  of 
120  VAC;  each  is  capable  of  supplying  750  W.  Initially,  the  voltage 
supply  was  supervised  manually  using  a  voltage  regulator.  Thermal 
control  was  good,  but  not  adequate  for  precise  measurements. 
Consequently,  we  are  currently  evaluating  a  controller  (Cole  Parmer 
2 18 6-2 OA)  which  responds  to  feedback  from  a  thermistor;  this  will 
require  one  more  probe  in  the  experimental  chamber,  but  potential 
improvements  in  the  control  of  the  chamber  temperature  warrant  this 
change.  Similar  control  is  expected  for  a  fluid-filled  chamber. 

Flash  system 

In  essence,  a  flash  system  consists  of  a  source  of  radiant 
energy  such  as  a  lamp,  possibly  an  aperture  plate,  a  specimen,  and 
a  temperature  measurement  system  on  the  non-illuminated  side  of  the 
sample.  Here  we  discuss  the  energy  source.  We  selected  an  EG&G 
linear  xenon  flash  tube  (0.4  cm  diameter  and  11.5  cm  long)  which  is 
rated  to  deliver  a  240  J  explosion  energy  at  a  pulse  duration  of  < 
135  /xs.  The  custom  circuit  for  the  lamp  source  is  shown  in  Figure 
15;  note  that  the  136  fiF  capacitance  is  actually  built  up  using  a 
series  of  smaller  capacitors.  The  calculated  energy  output  from 
this  circuit  is  80  J  based  on  a  1000  Amp  input  current.  The 
operation  of  the  lamp  is  straightforward.  In  the  non-ionized  state, 
the  flash  tube  has  a  high  impedance,  which  causes  the  current  from 
the  power  supply  to  flow  into  the  charging  capacitor.  The  instant 
the  trigger  is  released,  the  voltage  across  the  tube  increases  and 
xenon  atoms  within  the  tube  ionize.  This  causes  a  rapid  drop  in  the 
impedance  of  the  tube.  In  a  very  short  period,  enough  atoms  are 
ionized  to  allow  the  current  to  flow  from  the  capacitor  through  the 
flash  tube.  As  this  continues,  an  arc  expands  to  fill  the  tube. 
However,  in  a  short  period  most  of  the  current  from  the  capacitor 
is  released,  and  the  atoms  deionize  and  stop  conducting.  This 
allows  the  capacitor  to  recharge. 

To  better  focus  the  energy  from  the  lamp  (which  is  emitted 
radially)  onto  the  specimen,  a  linear  cylindrical  mirror  is  placed 
behind  the  lamp  (Figure  14)  and  the  surface  of  the  specimen  is 
darkened  with  carbon  black  or  India  ink.  As  mentioned  above,  the 
specimen  is  also  illuminated  through  a  square  aperture  (recall 
Figure  8a)  in  order  to  heat  only  a  portion  of  the  surface  —  this 
allows  the  in-plane  components  of  the  thermal  diffusivity  to  be 
determined.  The  aperture  will  be  mounted  within  the  device  in  order 
to  place  it  as  close  to  the  specimen  surface  as  possible.  Thus, 


although  the  total  energy  from  the  lamp  can  be  calculated,  that 
portion  which  reaches  the  sample  must  be  determined  by  placing  an 
instrumented  (with  thermocouples)  blackened,  high  conductivity 
metallic  plate  of  known  mass  and  specific  heat  just  in  front  of  the 
specimen  and  discharging  the  lamp;  the  high  conductivity  ensures 
nearly  uniform  heating.  The  energy  absorbed  is  found  via  0=Ae  = 
CpfliAr  (recall  equation  25)  . 

Temperature  measurement 

Two  sets  of  temperatures  must  be  measured  on-line;  the  ambient 
temperature  in  the  central  region  of  the  chamber,  and  the 
temperature  at  three  orthogonal  locations  on  the  non-illuminated 
surface  of  the  specimen  (Figure  8)  .  For  reasons  of  cost  effective¬ 
ness  and  general  utility,  thermocouples  are  currently  being  used 
(below  we  document  design  considerations  based  on  possible  IR 
camera  measurements) .  The  ambient  temperature  is  easily  measured 
with  a  T-type  thermocouple.  In  contrast  to  the  pilot  studies 
mentioned  above,  the  surface  temperature  is  measured  at  three 
orthogonal  locations  (Figure  8)  using  E-Type  thermocouples  that 
have  response  times  of  <  0.05  s,  a  nominal  bead  size  of  the  order 
of  0.033  cm,  and  a  wire  diameter  of  about  0.0254  cm.  These  three 
thermocouples  are  each  placed  in  a  0.035  diameter,  2.54  cm  long 
glass  tube  (to  provide  stability,  insulation,  and  electrical 
isolation)  ,  which  in  turn  is  mounted  in  a  1.905  cm  Delrin  disk. 
Mounted  in  this  way,  the  thermocouples  form  an  equilateral  triangle 
of  4.5  mm  per  side;  this  corresponds  to  a  length  slightly  larger 
than  the  size  of  the  current  aperture  through  which  the  specimen  is 
illuminated.  The  center  thermocouple  is  placed  at  the  center  of  the 
specimen,  which  coincides  with  the  center  of  the  illumination 
region,  in  order  to  measure  the  out-of-plane  diffusivity  whereas 
the  other  two  measure  the  in-plane  diffusivities.  The  Delrin  disk 
is  mounted  on  a  moveable  aluminum  bracket,  which  allows  the  three 
thermocouple  beads  to  be  simultaneous  brought  into  contact  with  the 
surface  (with  the  appropriate  alignment) .  The  placement  of  the 
probe  is  accomplished  via  a  micrometer.  Because  the  surface 
temperatures  will  be  measured  at  multiple  states  of  finite  strain 
(each  at  mechanical  equilibrium) ,  there  is  a  need  to  re-adjust  the 
position  of  the  probe  for  each  strain  state.  Moreover,  because  the 
temperature  measurement  is  a  function  of  the  contact  force  ( i . e . , 
contact  conductance) ,  a  calibrated  spring  will  be  mounted  on  the 
probe  to  allow  repositioning  at  a  fixed,  prescribed  indentation 
load. 

Although  not  currently  used,  we  evaluated  the  potential  advan¬ 
tages  of  infrared  temperature  measurement,  that  is  IR  cameras.  IR 
cameras  measure  entire  temperature  fields  r(X,t) ,  and  thus  may  be 
used  to  estimate  the  specific  heat  (recall  equation  25)  as  well  as 
provide  data  needed  for  a  Marquardt-based  determination  of  the 
components  of  (equation  27) .  The  decision  to  use  IR  cameras  is 
not  a  trivial  one,  however,  as  we  shall  show. 

It  is  well  known  that  all  matter  warmer  than  0  K  emits 
electromagnetic  radiation.  An  ideal  emitter,  or  black-body,  will 
have  a  spectral  emission  described  by  Planck's  distribution,  with 
a  total  emissive  power  proportional  to  .  Thus,  one  can  determine 
the  temperature  of  an  object  by  measuring  the  power  emitted  from 


it.  Wien's  Displacement  Law  further  states  that  the  wavelength  of 
maximum  spectral  emission  is  related  to  temperature  by  A  T  - 
2897.6  Km,  where  is  the  wavelength  at  the  maximum  emissive 
power  and  T  is  the  absolute  temperature  of  the  surface.  For 
example,  the  temperature  during  heat-treatment  of  metals  is  often 
estimated  by  the  specimen  color:  for  example,  A^3,j=9.9,  7.8  and  0.5 
^m  for  T=293.2,  373.2,  and  5800  K,  respectively.  Disadvantages, 
however,  are  that  IR  temperature  measurement  is  expensive, 
inaccuracies  arise  if  the  emissivity  of  the  surface  is  unknown,  and 
the  response  time  for  full  field  interrogation  needs  improvement. 
Furthermore,  measurements  cannot  be  made  through  a  medium  which 
absorbs  IR  radiation,  such  as  water.  Yet,  IR  cameras  remains  an 
active  area  of  R&D  —  size,  minimum  discernible  temperature  (MDT) , 
and  scanning  rate  are  constantly  being  improved  —  and  the  desire 
to  measure  Cp{C,T)  suggests  that  they  must  be  considered  as  a 
candidate  system.  Table  1  contains  information  on  the  commercially 
available  cameras  considered  in  the  summer  of  1995.  Of  these,  the 
Cincinnati  Electronics  and  AGEMA  models  were  examined  in  our 
laboratory  using  the  prototype  load-frame  and  flash  system  (recall 
Figure  4) .  Results  suggested  that  before  these  cameras  could  be 
used  to  their  full  potential  (a)  optional  lenses  would  be  needed, 
and  (b)  a  method  of  accurately  digitizing  the  full-field  signal 
would  need  to  be  developed.  The  optional  lenses  available  for  long 
distance  imaging  are  typically  between  $20  and  30k;  IR  lens  are 
much  more  expensive  than  standard  lens  since  they  are  made  in  much 
smaller  quantities.  Video  boards  with  enough  speed  and  buffer 
memory  to  record  data  in  real  time  and  with  adequate  resolution 
were  about  $  5k.  Hence,  it  appears  that  independent  of  unresolved 
technical  disadvantages,  an  entry  level  IR  system  would  cost  >  $ 
90k,  which  was  well  above  our  equipment  budget. 

Consequently,  we  then  evaluated  the  possibility  of  designing 
a  custom  temperature  measurement  system  using  discrete  IR  detectors 
(as  compared  to  the  128  x  128  or  256  x  256  arrays  of  typical  IR 
cameras)  .  That  is,  three  IR  detectors  could  be  used  in  an  array, 
and  thereby  provide  information  similar  to  that  given  by  the 
thermocouple  probe.  Given  the  expected  temperature  range  in  our 
studies  (and  hence  the  peak  emissive  power  wavelengths) ,  we  found 
that  Mercury-Cadmium-Telluride  (MCT)  and  Mercury-Manganese- 
Telluride  (MMT)  sensors  have  the  highest  D*  ("10^^cmHz°*®W^)  ,  a 
scaled  measure  of  the  IR  sensitivity  of  the  detectors.  Vendors  that 
were  contacted  were  Electro-Optical  Systems  (Phoenixville,  PA) , 
Brimrose  Corporation  of  America  (Baltimore,  MD) ,  and  EG&G  Judson 
Optoelectronics.  Brimrose  (MMT)  and  Electro-Optical  Systems  (MCT) 
were  willing  to  develop  unique  detectors  with  the  characteristic 
dimensions  needed  for  the  measurements,  average  costs  of  which  were 
quoted  to  be  about  $10  k.  EG&G  was  not  willing  to  develop  unique 
sensors  at  a  competitive  cost.  The  draw-backs  of  the  discrete 
sensors,  however,  are  (a)  the  need  to  design  and  develop  unique  IR 
optics,  (b)  the  lack  of  a  priori  knowledge  of  their  thermal 
sensitivity  and  calibration,  and  (c)  cost.  Electro-Optical  Systems, 
for  example,  suggested  purchasing  a  single  sensor  at  $2k  to  test 
its  applicability.  This,  in  turn,  would  have  required  purchasing  IR 
optics,  estimated  at  $5k  (including  optical  mounts).  A  >  $7  k 
expenditure  to  evaluate  a  sensor  was  not  justifiable. 

Nonetheless,  a  design  of  the  requisite  optics  was  developed. 


As  with  visible  optics,  there  are  both  transmission  and  reflective 
IR  optics.  Transmission  IR  optics  follow  the  same  principles  as 
optics  used  in  the  visible  spectrum,  but  require  exotic  materials 
such  as  germanium  (Ge)  and  zinc  selanide  (ZnSe) .  These  materials 
have  physical  properties  that  make  them  generally  difficult  to 
machine  and  use;  Ge,  for  example,  is  opaque  to  visible  light,  which 
makes  alignment  difficult.  An  advantage  in  fabrication,  however,  is 
the  long  wavelength  in  the  IR  spectrum  —  optical  tolerances  are 
often  quoted  in  terms  of  fractions  of  a  wavelength,  usually  the 
wavelength  in  the  center  of  the  spectrum  of  interest;  for  example, 
a  quarter  wavelength  tolerance  for  a  HeNe  laser  at  632  nm  requires 
more  precise  machining  than  for  an  IR  component  at  10  lira.  Optical 
design^  was  led  by  Dr.  Henry  Cox,  a  retired  optical  designer  from 
Annapolis  MD,  who  donated  his  time  to  the  project.  His  proposed 
design  (Figure  16)  consists  of  two  parabolic  mirrors,  with  focal 
lengths  of  240  and  300  mm,  respectively.  The  mirrors  were  to  be 
machined  from  6061-T6  aluminum  alloy,  and  the  optical  faces  gold 
plated.  Accurate  manufacture  of  the  optical  elements  is  critical. 
Quotes  for  the  mirrors  were  obtained  from  Janos  Technology  ($3  k) 
and  the  Center  for  Applied  Optics  at  the  University  of  Alabama  at 
Huntsville  ($  2.8  k) .  Material  to  support  the  mirrors  would  be 
machined  at  UMBC. 

In  summary,  the  expense,  uncertainty,  and  need  to  measure 
temperatures  on  submerged  specimens  led  to  the  use  of  E-type 
thermocouples.  Nonetheless,  the  potential  advantages  of  having  the 
entire  temperature  field  on  the  non-illuminated  face  remains 
attractive,  and  IR  measurement  techniques  will  continue  to  be 
evaluated  as  the  technology  improves. 

Data  Acquisition  and  Control 

The  description  of  the  primary  sub-systems  and  theoretical 
framework  highlights  the  need  for  adequate  acquisition  of  the 
required  data.  Depending  on  the  particular  test  protocol  employed, 
as  many  as  6  independent  channels  of  information  will  need  to  be 
recorded  via  the  A/D  board  (two  forces,  and  four  temperatures)  and 
4  surface  markers  will  have  to  be  tracked  using  the  video  system 
(i.e.,  4  (X2^,X2)  pairs).  Of  course,  the  strain  and  diffusivity 
related  temperatures  need  not  be  measured  simultaneously.  Similar 
to  the  signal  conditioner/amplif iers  for  the  force  transducers  (SA- 
4  card) ,  separate  amplifiers  (Omega  OmniAmp  III)  were  purchased  to 
amplify  the  thermocouple  signals  1000  fold.  Amplification  improved 
the  signal  to  noise  ratio  considerably.  Specifically,  the  load  and 
temperature  signals  are  read,  via  a  DT707  external  jumper  board,  by 
a  standard  12  bit  A/D  board  (DT  2837) ;  this  board  is  specifically 
designed  for  acquisition  of  thermocouple  data  and  thus  is  equipped 
with  a  cold  junction  compensation  circuit  at  the  0  channel.  The 
theraocouples  occupy  channels  1-4,  and  the  two  load  cells  channels 
positions  5  &  6,  respectively.  The  computer  platform  is  a  90  MHz 
Compaq  Pentium  personal  computer,  and  a  custom  program  allows  one 
to  vary  the  sampling  rate  and  gain  per  channel  independently. 


^  See  Modern  Optical  Engineering  by  W.  Smith  (McGraw-Hill, 
1990)  . 
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'ttlG  COHipU'td^XOndX  tlXlUG  fox'  O^Oh  X^OlTclf XOn  SlXsO  XncX'GclSGS 

For  exampXe  a  modeX  using  21  x  21  x  9  nodes  takes  approximately  20 
seconds  per  iteration  on  a  Sun  SparcStationS,  whereas  a  61  x  61  x 
9  node  model  takes  20  minutes  per  iteration.  Thus  when  the  number 
of  iterations  are  as  high  as  50  (for  noisy  data)  it  takes 
significant  time  to  reduce  the  data.  ConseguentXy,  we  are  currently 
developing  models  that  have  a  higher  concentration  of  nodes  where 
needed  (i.e.,  in  regions  of  steep  temperature  gradients),  but  fewer 
nodes  overall.  The  governing  equations  must  be  transformed  so  that 
we  may  relate  the  nodal  position  as  a  given  function. 

CONCLUSIONS 

Originally  we  proposed  to  use  an  existing  experimental  device 
to  perform  "standard"  tension-torsion  thermoelastic  tests  (e.g., 
see  Treloar,  1975)  on  a  number  of  elastomers  and  fiber-reinforced 
elastomeric  composites.  That  is,  the  emphasis  was  proposed  to  be  on 
data  collection  and  constitutive  formulations  for  a  specific  set  of 
materials.  Based  on  early  discussions  with  ARL-MO  personnel, 
however,  it  became  clear  that  the  focus  should  be  changed  to  the 
design  of  a  new  experimental  device  for  performing  tests  on  small 
samples  of  candidate  materials  as  they  become  available.  Hence,  we 
have  designed  a  theoretically-motivated  multiaxial  test  system  that 
is  capable  of  collecting  thermoelastic  and  thermophysical  data  on 
small,  thin  sheets  of  elastomers  that  are  subjected  to  a  variety  of 
finite  strains,  temperature  levels,  and  environmental  conditions. 
We  submit  that  this  device  represents  a  significant  advance  in 
elastomer  testing,  and  it  will  allow  cutting-edge  thermomechanical 
constitutive  relations  to  be  rigorously  formulated.  Many  different 
design  options  were  evaluated,  theoretical  concepts  were  developed, 
and  numerical  simulations  were  used  to  identify  preferred  tests. 
The  device  and  much  of  the  associated  control  and  data  reduction 
software  is  now  complete. 
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Table  1.  IR  cameras 


company 

model 

sensor 

cost 

($k) 

MDT  array  size  scan  rate 

(K)  (pixels)  (Hz) 

AMBER 

Radiance  1 

InSb 

59 

<0.025 

256x256 

60 

C-E 

IRRIS 

InSb 

80 

<0.025 

256x256 

50-150 

IRC 

InSb 

30 

<0.025 

160X120 

50 

I-M 

700  series 

HgCdTe 

60 

<0.2§30Hz 

256x200 

7866 

AGEMA 

450 

HgCdTe 

27 

0.1°C 

100X140 

30 

470 

HgCdTe 

80 

0.1°C 

100x140 

30 

900 

HgCdTe 

140 

<0.1°C 

200X136 

30 

FSI 

Prism 

PtSi 

32 

0.1°ce30°C 

320X244 

60 

816 

HgCdTe 

80 

0.05 

348x240 

30 

ISI 

VideoTherm 

vidicon 

21 

0.15®ce25®C  272  lines 

30 

C-E  is  Cincinnati  Electronics 

I-M  is  Infra-Metrics 

FSI  is  FLIR  Systems,  Inc. 

ISI  is  ISI  Group,  INC. 


Figure  1.  Sctiemetiic  diegram  of  'the  decomposition  of  a  thermoelastic 
deformation  into  that  part  due  to  heating  alone  (F-,)  and  that  part 
due  to  isothermal  stretching  (1^2) .  ^  ^ 


Figure  2  •  Schema  of  the  uniform  flash  method  proposed  by  Parker  et 
measuring  the  out-of-plane  diffusivity  (panel  a)  and  the 
modification  of  Malliet  et  al  for  measuring  2  components  of  a  (b)  . 
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Figure  3 .  Simulations  showing  potential  temperature  histories  o, 
the  non-illuminated  face  of  stretched  NGR  sample. 


Prel;minary  biaxial  load  frame  constructed  to  perform 
feasibility  studies  in  the  Spring  of  1995.  pertorm 
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Figure  7.  Finite  element  calculations  of  non-uniform  stress  (panel 
a)  and  strain  (panel  b)  distributions  that  are  induced  by  5 
individual  point  loads  per  side  (equibiaxial  stretch  of  25%) .  Note 
the  homogeneity  over  the  central  region. 
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Figure  8.  Schema  of  our  proposed  flash  method  for  measuring  the  3 
diagonal  components  of  the  thermal  diffusivity:  panel  a  shows  the 
set-up  and  panel  b  the  placement  of  the  thermocouples  relative  to 
the  markers  used  for  measuring  the  biaxial  finite  strain. 
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Time  vs.  Temp,  on  Rear  Face 


phi=((t-tamb)*VCp*rho)/q 
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Figiure  11.  3-d  heat  transfer  simulation  results 


Figure  I3a.  Oblique  view  of  the  load  frame  (LF)  and  environmental 
chamber  (EC) . 
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Figure  15.  Circuit  diagram  for  the  flash  tube. 
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Figure  17.  Dependence  of  the  SDecifio 
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ABSTRACT 

Investigations  comprising  of  experimental  observations  and  analytical  study  have  been 
systematically  executed  to  evaluate  a  particular  composite  system  which  has  the  potential  for 
structural  armor  applications.  The  material  being  a  thick  laminate,  is  a  random  mode  stacking  of  bi¬ 
axial,  plain  weave  S2-Glass/S2-Glass  fiber  yams  with  polyester  resin  (Cycom  4102)  as  binder 
material,  thirty  two  percent  by  weight.  Attention  is  focused  to  understand  the  three  dimensional 
macroscopic  behavior  of  the  material  under  quasi-static  and  dynamic  compressive  loading 
conditions.  Quasi-static(low  to  medium  strain  rates)  experimental  investigations  were  carried  out  on 
an  MTS  test  system  while  the  dynamic  response  was  observed  using  the  Split  Hopkinson  Pressure 
Bar  test  equipment.  The  uniformity  of  the  material's  response  is  verified  for  different  sized  and 
shaped  specimens.  The  material  has  been  found  to  be  transversely  isotropic  and  its  three 
dimensional  stiffness  matrix  has  been  determined.  The  material's  strength  and  strain  rate  effects  in 
the  region  of  1  O'*  -  10^  sec'*  are  also  presented.  The  materials  failure  stress  has  been  found  to  be 
higher  under  dynamic  loads.  The  failure  mechanism  for  the  material  in  direction  perpendicular  to 
ply  layup  is  fairly  complex  as  compared  to  the  other  two  orthotropic  directions.  Effort  has  been 
made  to  characterize  the  damage  modes  by  investigating  fractured  surfaces. 


INTRODUCTION 


In  continuing  search  for  light  weight  materials  of  desirable  strength  and  stiffness, 
considerable  effort  has  been  made  in  the  past  six  decades  in  the  technological  development  of  fiber- 
reinforced  materials‘'\  Their  low  specific  gravities,  the  high  strength-weight  and  modulus-weight 
ratios  of  these  composite  materials  are  markedly  superior  to  those  of  metais^‘\  A  potential  for  weight 
savings  with  fiber  reinforced  composites  exists  in  many  engineering  fields,  implementing  them  for 
actual  use  would  necessitate  rigorous  design  practices  and  viable  process  developments  based  on  the 
understanding  of  their  mechanical  and  physical  characterization.  Composite  structural  elements  are 
now  in  a  variety  of  components  for  automotive,  aerospace,  marine  and  architectural  structures  in 
addition  to  consumer  products  such  as  skis,  golf  clubs,  tennis  rackets,  etc.  One  of  the  earliest  and 
ever  growing  application  area  is  in  the  military.  It  is  well  known  that  performance  and 
maneuverability  of  vehicles  (airborne  or  ground)  depend  substantially  on  weight.  High  specific 
strength  and  high  specific  stiffness  characteristics  of  the  composites  encouraged  the  technical 
community  to  endeavor  for  its  use  in  protective  armor. 

The  first  use  of  Glass  fiber-reinforced  plastic  for  armor  applications'*’  occurred  with  the 
development  of  Doron“  which  was  employed  by  the  US  Marines  during  the  Korean  conflict.  With 
the  development  of  Kevlar'’,  its  applications  ranged  from  Personal  Armor  systems  such  as 
fragmentation  vest  &  helmets  to  spall  liners  within  the  Ml  13  APC  and  Bradley  Fighting  Vehicles. 
S2-Glass  reinforced  plastic  composite  is  now  under  consideration  for  combat  vehicle  hulls.  For  such 
a  load  bearing  armor  application  it  is  but  essential  to  evaluate  the  structural  integrity  and  dynamic 
response  of  these  composites.  It  is  imperative  that  a  complete  understanding  of  the  mechanical 
behavior  of  the  composite  system  is  made  known  and  its  behavior  under  different  types  of  loading 
condition  be  predicted  through  simulations,  before  the  material  could  be  implemented  for  its  desired 
application.  The  application  is  such  that  it  is  also  necessary  to  understand  the  residual  strength  and 
performance  of  the  composite  system  upon  ballistic  impacts. 

Chou  and  Ko  have  illustrated  that  woven  composites  have  proven  to  provide  better 
stability  in  terms  of  strength  and  stiffness  and  extensibility  than  that  of  a  unidirectional  prepreg.  This 
has  interested  researchers  like  Chen  &  Kuo^**’,  Ishikawa'^"*’,  Hashin,  Rosen,  Naik'^'^’et  al  to  conduct 
studies  on  2-D  woven  SiC/SiC  and  C/SiC  composites  and  have  formulated  their  effective  stiffness 
and  strength.  The  above  works  are  focussed  on  analytical  and  numerical  results  with  plane  stress 
and  plane  strain  conditions.  The  strain  in  the  thickness  direction  observed  in  such  composites  are 
significant  and  cannot  be  neglected  from  the  stress-strain  relation. 


®A  laminate  of  fiber  E-glass  and  polyester  resin. 
'’A  high  tenacity  fiber  developed  by  Dupont. 


Naik  and  Ganesh'"*'  carried  out  some  shear  strength  studies  on  E-glass/epoxy  resin  plain 
weave  fabric  laminates.  The  experiments  were  confined  to  tensile  loading  at  slow  rates.  Karayaka 
and  Kurath'"’  conducted  failure  studies  on  Graphite/Epoxy  plain  weave  laminate  under  compression, 
tension  and  three  point  flexural  loading.  Kumar  and  Garg  investigated  failure  modes  in 
unidirectional  E-glass-epoxy  cross  ply  laminates  under  dynamic  compression  and  characterized  them 
fractographically.  Again  the  samples  investigated  in  the  above  cases  were  from  laminates  whose 
thickness  were  less  than  3mm.  Dandekar,  Green  and  Beaulieu^^’  had  carried  out  studies  on  thick 
S2-GRP  to  determine  the  elastic  constants  through  ultrasonic  techniques  and  had  characterized  the 
material  to  be  transversely  isotropic.  They  also  performed  quasi-static  tests  in  compression  &  tension 
along  with  some  high  strain  rate  experiments  in  the  thickness  direction. 

As  mentioned  previously,  identifying  the  onset  of  damage  and  its  evolution  while  the 
composite  is  under  load  or  upon  impact  is  being  further  evaluated.  Understanding  the  damage 
mechanics  is  very  useful  as  it  can  aid  to  predict  the  residual  strength  of  the  composite  and  its 
response  to  subsequent  loads/impacts.  Damage  development  studies  in  unidirectional  graphite/epoxy 
cross-ply  laminates  were  carried  out  by  Daniel  and  Lee“^\  Damage  was  characterized  by  the 
classical  decreasing  stiffness  notation  and  its  correlation  to  fiber  failure  and  matrix  cracking.  Shan, 
Pulvinage^*"*’  et  al  reported  analysis  and  modelling  of  damage  behavior  of  two-dimensional  woven 
SiC/SiC  composites.  Their  work  was  supported  by  numerical  results.  Degradation  of  the  elastic 
properties  of  matrix,  fiber  and  interfacial  bonding  was  considered  as  a  measurement  for  the  evolution 
of  damage.  This  classical  approach  of  relative  variation  of  elastic  properties  as  an  indicator  of 
damage  cannot  predict  the  origin  of  damage  and  its  location.  Acoustic  emission  detection  technique, 
upon  careful  monitoring  can  provide  us  with  information  pertaining  to  the  origin,  evolution  and 
location  of  damage.  Sachse  and  Kim^‘^\  and  Liao  and  Targ^‘*\  Gardiner  and  Pearson^’’^  have  all 
contributed  to  use  of  acoustic  emission  in  understanding  material  properties  variation  and  monitoring 
of  damage  occurring  under  static  loading.  Their  work  was  focussed  on  graphite/epoxy  unidirectional 
laminates  and  involved  the  traditional  acoustic  emission  parameters  such  as  amplitude,  count, 
energy,  time,etc.,  to  quantify  the  failure  mechanisms.  Conventional  AE  methods  employed  resonant 
sensors  and  narrow  band  pass  filters.  This  type  of  measurement  discards  much  of  the  information 
available.  With  recent  advances  in  technology,  it  is  now  possible  to  perform  wave  based  AE 
measurements,  along  with  high  fidelity  and  broad  band  sensors.  This  allows  to  analyze  the  damage 
occurrences  and  its  evolutions  based  on  the  physical  dynamics  of  the  micro-structure  of  the 
composite  material.  Gorman^'*’  in  his  paper  presents  a  comparative  illustration  between  the  old  AE 
method  and  the  new  wave  based  AE  method.  The  damage  evaluation  experiments  using  the  wave 
based  acoustic  emission  equipment  is  presently  being  carried  out  in  the  laboratory. 


EXPERIMENTAL  PROCEDURE 


(i)  Material 

The  material  constituents  in  the  composite  are  S2-Glass  in  the  fiber  form  and  Polyester 
resin(CYCOM  4102)  as  the  binder  or  matrix.  The  construction  of  the  composite  is  as  follows: 
Yam :  The  yam(see  Fig.  1)  is  a  collection  of  parallel  continuous  ends  of  fibers  of  average  diameter 
of  10  microns.  The  fibers  are  so  gathered  to  give  the  yam  an  approximate  elliptical  cross-section  with 
the  ratio  of  major  to  minor  axis  being  around  10  .  Yam  typically  consists  of  approximately  10‘° 
fibers  bundled  together.  A  binder  material  is  distributed  in  the  inter  fiber  space  to  the  desired  volume 
fraction. 

Lamina :  The  lamina(see  Fig.2)  is  a  two  dimensional  plain  weave  balanced  constmction  of  the  fiber 
yams  in  the  two  mutually  perpendicular  direction.  'Wrap  yam'  is  termed  for  the  yarn  in  one  weave 
direction  and  'fill  yam'  for  the  yam  in  the  other  perpendicular  weave  direction.  The  woven 
constmction  is  then  impregnated  with  the  polyester  resin  in  the  inter-yam  space  so  as  to  obtain  a 
lamina  of  uniform  thickness.  The  lamina  measured  O.bSmm  in  thickness. 

Laminate  :  The  laminate(see  Fig.3)  is  a  random  phase  mode  multi-layer(42  ply)  stacking  of  such 
laminas  such  that  the  wrap  and  fill  yams  align  in  the  same  directions.  The  stacked  lamina  is  then 
cured  under  desired  pressure  and  temperature  cycle  in  an  autoclave.  The  material  was  in  plate 
form(5lcmx5lcm)  having  42  plies  and  measured  26.4mm  in  thickness.  The  materials  density  was 
found  to  be  1.92  ±  0.03  gm/cm\  The  resin  content  was  32  ±  2  percent  by  weight.  The  properties 
of  CYCOM  4102  along  with  the  fiber  volume  fraction  in  the  material  construction  are  provided  in 
Table  1. 


Fig.l:Yam  Fig.2:  Lamina  Fig.3  :Laminate 

The  fiber  orientation  and  ply  lay-up  suggests  that  the  material  mechanical  response  would  be 
orthotropic.  Three  orthotropic  direcdons(see  Fig.4)  were  identified  namely,  [100]  Direction  -  along 


the  axis  of  the  wrap  yam,  [OlOJ  Direction  -  along  the  axis  of  the  fill  yam  and  [001]  Direction- 
perpendicular  to  the  lamina  plane 


2 


The  composite  was  observed  not  to  absorb  moisture.  In  machining  the  test  specimens  Kennametal 
KC730  Titanium  Nitride  and  Diamond  Carbide  cutting  tools  were  employed. 


Table  1 


Density 

gm/cm^ 

Modulus 

GPa 

Strength 

MPa 

Weight 

Percent 

Volume 

Percent 

S2-Glass  Fiber 
(10-12  microns  dia.) 

2.49 

85.5 

4585 

68 

52.37 

Polyester  Resin 
(CYCOM4102) 

1.3 

1 

3.6 

78.5 

32 

47.63 

(ii)QUASI-STATIC  COMPRESSION  TESTS 

The  material's  mechanical  response  in  compression  was  investigated.  The  quasi-static(  low 
to  medium  strain  rates)  compression  experiments  were  carried  out  in  the  MTS  material  test 
equipment.  The  MTS 809  is  a  servo-hydraulic  controlled  axial/torsion  test  equipment  with  a 
250kN/2200Nm  static  load/torque  capability.  The  controller  to  the  MTS  is  interfaced  with  a 
computer  and  Teststar  SX-790.10  is  the  software  module  used  to  program  the  compression 
experiments.  The  hydraulic  actuators  of  the  MTS  can  be  programmed  to  move  linearly  in  a 


continuous  range  from  lO"^  to  3  in/sec.  Force  measurements  were  made  by  means  of  load  cell  on  the 
I^rrS  and  the  strains  were  monitored  using  strain  gages  and  the  MTS  LVDT.  The  strain  gage  data 
was  used  in  the  strain  region  of  interest.  In  higher  strain  regions  the  LVDT  data  along  with  the 
correction  factor  was  used.  The  force,  displacement  and  the  strain  readings  were  acquired 
continuously  through  the  data  acquisition  built  in  the  controller  computer .  It  had  the  capability  of 
collecting  data  every  0.001  second.  The  load  cell  and  the  LVDT  were  calibrated  and  the  error  in  the 
readouts  were  less  than  2  percent  in  the  range  of  200kN  force  and  25  mm  displacement. 

Preliminary  experiments  were  conducted  to  investigate  the  shape  and  size  effects  on 
specimen  as  the  material  was  not  microscopically  homogeneous.  Two  types  of  geometries  were 
selected  i.e.circular  cross-section  and  square  cross-section.  Two  different  sizes  were  chosen  for  each 
geometry.  Fig.  4  gives  the  details  of  the  specimens  tested  for  this  purpose.  In  all  of  the  samples , 
three  types  of  specimen  orientation  were  machined  such  that  the  loading  axis  corresponded  to  the 
three  orthotropic  directions.  Table  2  summarizes  the  results  and  Figs.5  &  6  illustrate  the  material's 
directional  elastic  responses.  All  specimens  were  loaded  in  compression  under  a  constant  strain  rate 
of  10"*  sec  ‘.  Teflon  sheet(0. 127mm  thickness)  was  used  at  the  specimen  ends  for  lubrication. 

In  compression  tests  to  evaluate  the  elastic  constants'  of  the  material ,  specimens  with  square  cross- 
section  measuring  12.7x12.7  mm  and  25.4  mm  long,  were  chosen  .  In  addition  to  the  three 
specimen  orientation  chosen  earlier ,the  [110]  oreintation  specimen  were  also  used.  Two  element 
rosettes(KYOWA  Electronic  Instruments  Co.,  Ltd.,  Japan)  of  5  mm  gage  length  were  mounted  on 
3  adjacent  vertical  surfaces.  The  strain  readings  from  the  four  gages  on  adjacent  faces  would  describe 
the  axial  strain  in  the  loading  (one  orthotropic)  direction  and  its  corresponding  transverse  strains  in 
the  other  two  orthotropic  direction.  The  measured  elastic  constants  from  these  test  are  given  in 
Table3. 


1.0"xl.0"x2.0'’ 

Fig  4.  Specimen  samples 


Next  the  material’s  responses  for  different  strain  rates  i.e.  lO'^sec*  ,  10'^  sec‘‘  and  10°  sec‘‘, 
were  investigated  in  all  three  orthotropic  directions.  Fig.7-  10  illustrate  the  results  of  these  tests. 

(iii)  DYNAMIC  RESPONSE  EXPERIMENTS 

The  dynamic  response  of  the  material  was  investigated  under  high  strain  rate  compressive 
loading  by  means  of  split  Hopkinson  Pressure  Bar  test  equipment.  The  impacting  bar  was  56  cms 
long  which  resulted  in  a  loading  pulse  of  about  220  micro-seconds.  The  impacting  bar  is  propelled 
out  of  a  gas  gun,  using  compressed  nitrogen  gas.  The  magnitude  of  the  loading  pulse  was  a  function 
of  the  firing  pressure.  Strain  rates  upto  10^  was  attained  in  these  experiments.  The  incident  and 
transmitted  bars  diameter  was  0.99  inch.  By  means  of  strain  gages  on  the  incident  and  transmitted 
bars,  the  strain-time  histories  of  the  three  pulses,  incident,  reflected  and  transmitted  were  recorded 
on  a  Nicolet  440  digital  oscilloscope  with  12  bit  digitization  and  sampling  rates  as  low  as  0.1  micro¬ 
second.  The  nominal  stress  and  strain  in  the  specimen  upto  any  time  duration  t,  was  determined  by 
using  standard  relations  given  below,  for  this  technique. 
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However,  information  within  the  strain  regions  of  2%  cannot  be  considered  accurate  in  this 
technique.  The  specimen  were  cylindrical  disks  and  two  different  sizes  were  tested  in  the  [001] 
directions.  One  specimen  size  was  measuring  19  nun  in  diameter  and  9.6mm  thick  and  the  other  was 
measuring  12.7mm  in  diameter  and  6.35mm  thick.  Specimen  samples  in  the  [100]  and  [010] 
orthotopic  directions  were  of  the  same  size  being  19mm  in  diameter  and  9.6mm  thick.  Table  4  gives 
the  details  of  the  specimens  and  the  average  strain  rates  applied  in  the  experiments. 


RESULTS  AND  DISCUSSION 
Quasi-Static  Response 

The  material's  response  did  not  significantly  change  due  to  size  or  shape  of  the  specimens 
tested,  although  there  was  about  12  percent  variation  between  maximum  and  minimum  moduli 
values  as  seen  in  Table  2.  This  variation  cannot  be  identified  as  a  size  effect  because  a  similar 
observation  is  made  in  Table  3  for  specimens  of  the  same  geometry  and  size.  This  variation  suggests 
that  the  nature  of  the  micro-inhomogeneities  within  the  material  which  is  inherent  due  to  the  manner 


in  which  the  material  was  fabricated.  The  materials  elastic  response  in  the  [100]  and  [010]  direction 
was  observed  to  be  very  similar  and  linear.  The  average  modulus  in  the  [100]  direction  was 
determined  as  29.67i3.65GPa  and  that  in  the  [010]  direction  was  determined  as  30.06±  1.74GPa. 
Dandekar^^\  et  al.  in  their  studies  had  reported  a  variation  in  the  modulus  from  26.3  to  47.5GPa  in 
[100]  and  from  24.5  to  43.0GPa  in 


Fig.5  :  Response  of  [100]  and  [010]  orientation  specimens 


Fig.6:  Response  of  [001]  oreintation  specimens 


Table  2 


SPECIMEN  DIMENSIONS  ELASTIC  RESPONSE  MODULI(GPa) 


iVITS_Xl 

QS^x  0.5”x  1.0” 

LINEAR 

29.62 

MTS_X2 

0,5”x  0.5”x  1.0” 

LINEAR 

29.26 

MTS_X3 

1.0”x  1.0”x  2.0” 

LINEAR 

29.48 

MTS.X4 

1.0”x  1.0”x  2.0” 

LINEAR 

31.11 

MTS^XS 

DIA.  0.5**  X  1.0*’ 

LINEAR 

33J3 

MTS_X6 

DIA.  0.5”  X  1.0” 

LINEAR 

31.42 

MTS.X7 

0.65”x  0.65”x  0.4” 

LINEAR 

26.67 

MTS.X8 

0.65*’x  0.65”x0.4” 

LINEAR 

28.77 

MTS_X9 

DIA.  0.75”x  0.4” 

LINEAR 

28.0 

MTS^XIO 

DIA.  0.75”x  0.4” 

LINEAR 

29.09 

MTS.Zl 

0.5**x  0,5”x  1.05” 

APP.  LINEAR 

13,05 

MTS^Z2 

0.5**x  0.5**x  1.05” 

APP.  LINEAR 

14.5 

MTS^Z3 

DIA.0.75  X  2.0” 

APP.  LINEAR 

11.25 

MTS.Z4 

DIA.0.75  X  2.0” 

APP.  LINEAR 

10,65 

MTS_Z5 

DIA.  0.5  *  X  1.05  * 

APP.  LINEAR 

11.76 

MTSJ6 

DIA.  0.5”  X  1.05** 

APP.  LINEAR 

12.73 

MTSJ7 

0.65”x  0.65”x  0.4” 

APP.  LINEAR 

16.11 

MTSJ8 

0.65’*x  0.65”x0.4” 

APP.  LINEAR 

11J8 

MTS^Yl 

0.5”x  0J”x  1.0” 

LINEAR 

28J2 

MTS^Y2 

0.5”x  0.5”x  1.0” 

LINEAR 

29.56 

MTS_Y3 

1.0**x  1.0”x  2.0” 

LINEAR 

30.12. 

MTS_Y4 

1.0”x  1.0”x  2.0** 

LINEAR 

30.81 

MTS.Y5 

DIA.  0.5”  X  1.0” 

LINEAR 

31.02 

MTS^Y6 

DIA.  0.5”  X  1.0” 

LINEAR 

30.14 

MTS_Y7 

0.65”x  0.65’*x  0.4** 

LINEAR 

29.45 

MTS^Y8 

0.65”x  0.65**x0.4” 

LINEAR 

30.13 

MTS,Y9 

DIA.  0.75**x  0.4” 

LINEAR 

29.67 

MTS_Y10 

DIA.  0.75”x  0.4” 

LLNEAR 

31.38 

[100]  Orientation 


[001]  Orientation 


[010]  Orientation 


SPECIMEN 

MODULUS 

Gpa 

POISSONS  RATIO 

Failure 

Strain,  % 

Failure 

Stress,  MPa 

MTS^Xl 

28.99 

V,j=  0.12 

V,3=  0.34 

1.02 

208.5 

MTS_X2 

30.85 

V,j=:  0.11 

V,j=  0-J7 

1.12 

2323 

MTS_X11 

31.27 

V„*0.14 

V,3=  0.35 

0.87 

187.4 

MTS_X12 

32.78 

V,J=  0.12 

V,j=0J4 

1.17 

245J 

MTS_X13 

30.67 

V„*0.12 

V,j=0.3S 

1.04 

2283 

MTS^Yl 

28J2 

Vj,=  0.12 

Vj3=0.34 

0.98 

215.6 

MTS^Y2 

29.56 

0.10 

Vu=0.32 

1.06 

236.8 

MTS^Yll 

30.12 

0.13 

Va=0J8 

1.20 

2523 

MTS^Y12 

32.02 

V„=0.12 

Va=0J5 

1.16 

248.3 

MTS_Y13 

31.23 

V,,=  0.14 

Va=0J4 

1.11 

243.5 

MTS^Zl 

14.36 

Vj,=  0.18 

V3j=:0.16 

6.83 

620.6 

MTS,Z2 

10.67 

V„=0.17 

V3J=  0.19 

7.43 

680.3 

MTS_Z9 

12.62 

V„=0.19 

V3J=  0.18 

7.56 

690.6 

MTS^ZIO 

11.95 

V3,=  0.19 

V3J=  0.18 

731 

611.6 

MTSJll 

11.22 

V3,=  0.19 

0.18 

832 

703.1 

MTS.XY.l 

1337 

— — 

MTS.XY.2 

14.05 

— 

— 

— 

— 

MTS.XYJ 

13.83 

— 

— 

— 

— 

MTS^XY_4 

12.78 

— — 

— 

— 

— 

MTSJXYJ5 

1333 

— 

— 

— 

— 

[010]directions.  It  was  noticed  that  the  response  was  linear  elastic  until  failure  in  both  of  these 
directions.  Failure  is  characterized  by  multiple  delaminations  and  loss  of  structural  integrity. 

The  average  failure  stress  for  [100]  specimens  was  found  to  be  220.36±  28.88  MPa  and  that  for 
the  [010]  specimens  was  239.3±  18.35MPa  as  compared*  to  212.9±23.2  and  210.3±46.2  Mpa, 
respectively.  The  average  failure  strains  recorded  were  1.04  and  1.10  percent  in  the  [100]  and 
[010]  orientation  as  again  compared®  to  0.76±0.18%  and  0.89±0.04%,  respectively.  The 
Poissons  ratio  was  higher  in  the  [001]  direction  with  a  average  value  of  0.35. 
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The  material’s  response  in  the  [001]  specimens  was  very  different  from  that  of  the  [  100]  and 
the  [010]  specimens.  The  modulus  was  much  lower  and  the  average  value  was  12.67GPa.  Slight 
non-linearity  is  observed  during  loading  and  the  unloading  occurs  on  a  different  stress-strain  loci  and 
is  more  non-linear.  Only  the  loading  stress-strain  region  is  used  to  compute  the  modulus  in  this 
direction.  This  difference  in  the  loading  and  unloading  stress-strain  region  is  due  to  the  polymer 
matrix  dominance  in  this  direction  of  deformation  and  the  polymer  matrix  is  visco-elastic  in  nature. 
This  also  pertains  to  the  lower  modulus  in  this  direction.  The  failure  stress  and  strains  have  been 
found  to  be  much  higher  than  that  in  the  [100]  and  [010]  directions.  A  lower  dispersion  in  the  failure 
stresses  and  strains  is  observed  and  the  average  values  for  the  failure  stress  being  660,8±41.24Mpa 
and  that  for  failure  strain  being  7.49±0.83  percent. 

The  failure  mode  in  this  direction  is  much  more  complicated.  Most  of  the  specimens  failed 
in  the  mid  region  and  across  one  or  two  shear  planes.  Fiber  yam  fracture,  matrix  cracking  and 
delamination  were  all  coupled  to  result  in  a  catastrophic  failure.  Blackened  regions  were  observed 
in  the  shear  area  indicating  possible  temperature  rise(Fig.l2).  Figures  10-13  show  the  photo 
micrographs  of  failed  surfaces  of  specimens. 

The  material's  strain  rate  dependence  in  the  [001]  direction  is  much  more  than  that  observed 
in  the  [100]  and  [010]  directions  for  the  strain  rates  region  10"^  to  10®.  This  is  observed  from  fig7- 
fig.9. 


Fig.7;  Effect  of  strain  rate  in  [100]  orientation  specimens 
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Fig.  10:  Photomicrograph(20X)  of  a  failed  [001]  specimen  under  quasi-static 
compression.  Failure  has  occured  across  two  adjacent  shear  planes 


Fig.  1 1:  Photoniicrograph(12X)  of  failed  [001]  specimen’s  surface,  viewing  from  top 
illustrating  fiber  yams  mpture  along  the  shear  plane. 


Table  4  gives  a  listing  of  specimens  details  with  average  strain  rates  applied  in  the  Split 
hopkinson  pressure  bar  tests  and  the  failure  stress,  if  failure  occured. 

Table  4 


Specimen 

No. 

Orientation 

Specimen  Dimensions 

Strain 

Rates 

Failure  stress 

MPa 

SHPB-22 

[001] 

Dia.  =  0.50",  Th.  =  0.25" 

313 

— 

SHPB-17 

[001] 

Dia.  =  0.50",  Th.=  0.25" 

434 

— 

SHPB-18 

[001] 

Dia.  =  0.50",  Th.  =  0.25" 

483 

— 

SHPB-27 

[001] 

Dia.  =  0.50"  ,Th.  =  0.25" 

506 

— 

SHPB-16 

[001] 

Dia.  =  0.50",  Th.  =  0.25" 

690 

910 

SHPB-15 

[001] 

Dia.  =  0.50",  Th.  =  0.25" 

1182 

924 

SHPB-28 

[001] 

Dia.  =  0.50",  Th.  =  0.25" 

1610 

970 

SHPB-5 

[001] 

Dia.  =  0.75",  Th.  =0.375" 

418 

SHPB-25 

[001] 

Dia.  =  0.750",  Th.  =  0.375" 

516 

SHPB-9 

[001] 

Dia.  =  0.750",  Th.  =  0.375" 

610 

— 

SHPB-26 

[001] 

Dia.  =  0.750",  Th.  =  0.375" 

639 

— 

SHPB-29 

[001] 

Dia.  =  0.750",  Th.  =  0.375" 

772 

SHPB-13 

[001] 

Dia.  =  0.750",  Th.  =  0.375" 

924 

SHPB-8 

[100] 

Dia.  =  0.750",  Th.  =  0.375" 

148 

185 

SHPB-20 

[100] 

Dia.  =  0.750",  Th.  =  0.375" 

320 

245 

SHPB-21 

[100]  1 

Dia.  =  0.750",  Th.  =  0.375" 

290 

203 

SHPB-24 

[010] 

Dia.  =  0.750",  Th.  =  0.375" 

543 

270 

SHPB-23 

[010] 

Dia.  =  0.750",  Th.  =  0.375" 

348 

190 

SHPB-19 

[010] 

Dia.  =  0.750",  Th.  =  0.375" 

143 

180 

Figs.  14  and  Fig.  15  shows  the  dynamic  stress-strain  responses  in  the  [001]  orientation  for  0.5"  dia. 
and  0.75"  dia  disk  specimens,  respectively. 
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As  mentioned  earlier,  SHPB  data  in  small  deformation  region(e.g.  below  1.5%)  is  not  reliable.  As 
it  was  observed  previously  in  quasi-static  experiments,  the  specimen  size  has  no  or  very  little  effect 
as  seen  from  these  two  figures.  Strain  rate  effect  is  prominently  observed  in  both  sizes  of 
specimens.  Another  interesting  observation  is  the  strength  of  the  material  in  the  [001]  orientation 
during  dynamic  loading  is  much  higher.  The  average  observed  failure  stress  was  noted  to  be  around 
930MPa  at  a  strain  rate  of  about  1 182/sec.  The  failure  strain  was  around  9.5%.  In  the  ARL®’  report, 
the  failure  stress  was  significantly  different  and  was  reported  to  be  at  370MPa  for  a  strain  rate  of 
3000/sec.  None  of  the  larger  diameter  specimens  failed  as  the  stresses  and  strains  had  not  reached 
the  failure  values.  Three  of  the  smaller  diameter  specimens  failed  at  stresses  of  910MPa,  924MPa 
and  970Mpa  but  at  the  same  strain  of  around  9.5%.  The  specimen  failure  mode  was  similar  to  that 
in  the  quasi-static  experiments,  resulting  in  two  parts  upon  complete  fracture.  In  one  case,  where 
the  strain  rate  was  greater  than  1600/sec,  the  specimen  failed  with  multiple  fragments. 

The  [100]  and  [010]  oreintation  specimens  tested  on  Split  Hopkinson  Pressure  Bar  equipment  did 
not  show  much  difference  in  strength  from  that  in  the  quasi-static  tests.  As  the  failure  strains  for 
these  orientations  are  in  the  range  of  1-1.5%,  the  stress-strain  plots  could  not  be  considered  very 
reliable.  The  failure  stress  observed  was  in  the  range  of  180-270MPa.  Although  the  failure  stresses 
were  quite  similar  to  that  in  quasi-static  loading,  the  failure  mode  was  interestingly  different.  The 
specimen  circular  cross-section  changed  to  an  ellipse  with  the  major  axis  in  the  direction  of  the  ply 
layup.  Upon  reaching  a  maximum  limit,  delaminations  began  to  occur  and  resulted  in  a  failure.  The 
number  of  delaminations  occuring  depended  upon  the  strain  rate  and  stress  pulse  magnitude. 


CONCLUSIONS 

Experiments  were  performed  over  a  wide  strain  range  in  specimens  with  three  orientation. 
The  material  behavior  was  found  to  be  complex.  In  the  two  fiber  directions,  the  response  was  almost 
elastic  until  failure  with  almost  no  strain  rate  sensitivity.  Failure  was  observed  to  be  mostly  due  to 
delaminations.  In  the  thickness  direction,  i.e.  perpendicular  to  fibers,  moderately  large  visco-elastic 
deformations  were  observed.  The  elastic  constants  and  failure  stresses  were  significantly  strain  rate 
dependent.  The  failure  mode  was  a  combination  of  yam  fracture,  matrix  cracking  and  delaminations 
across  shear  planes.  Black  regions  were  observed  on  failed  surfaces  suggesting  considerable 
temperature  rise  during  failure. 
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